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1.  INTRODUCTION 


1.1  Development  of  Gas  Turbine  Blade  Cooling 

Aerospace  turboengines  present  a  demanding  challenge  to  many  heat  transfer 
scientists  and  engineers.  Designers  in  the  field  are  always  seeking  the  best  design  to 
transform  the  chemical  energy  of  the  fuel  into  the  useful  work  of  propulsive  thrust  at 
maximum  effectiveness  (D.  G.  Shepherd,  1972,  J.  H.  Horlock,  1985).  Based  on  the 
principles  of  thermodynamics,  the  inlet  temperature  and  pressure  of  the  gas  in  aerospace 
turboengines  must  be  increased  in  order  to  improve  the  thermal  efficiency  of  the 
turboengines.  As  a  result,  advanced  gas  turbines  are  designed  to  work  at  a  inlet 
temperature  of  1500-1700  K,  and  the  turbine  blades  must  operate  at  very  hostile 
environments.  This  hostile  working  condition  is  often  in  conflict  with  the  material 
capability  of  withstanding  these  high  temperatures  (E.  R.  Norster,  1969). 

Turbine  blade  cooling  was  first  considered  in  a  German  design  in  1935.  It 
presented  a  very  tough  but  important  problem  to  heat  transfer  researchers  and  scientists 
due  to  the  high  temperature  and  rotation  of  turbine  blades  (Bathie,  1984,  Cohen,  1987). 

There  are  basically  three  approaches  considered  by  the  thermal  engineer 
worldwide  to  solve  these  problems.  The  first  approach  is  to  develop  a  new  surface 
coating  that  can  withstand  these  hostile  environments  and  has  a  very  low  heat 
conductivity  to  protect  the  blade  surface.  This  approach  is  very  practical  and  valuable. 
However,  the  progress  in  this  approach  is  generally  slow.  Because  this  new  surface 
coating  must  be  subject  to  the  impingement  of  high  temperature  gas  and  the  vibration  of 
turboengines,  the  new  surface  coating  is  difficult  to  attach  onto  the  blade  surface  firmly. 
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The  second  approach  is  the  air  cooling  technology  that  is  the  primary  method 
used  in  the  blade  cooling  so  far.  The  air  cooling  technology  is  a  traditional,  convenient, 
and  effective  blade  cooling  means.  It  is  subdivided  into  the  following  three  kinds: 

(1)  Convection  cooling,  which  is  an  early  blade  cooling  means.  Cooling  air  from 
the  compressor  flows  into  the  inner  channels  in  the  blade  from  the  blade  root,  and  flows 
out  from  the  blade  top.  The  spent  air  is  discharged  and  mixed  with  the  main  gas  flow  of 
the  next  stage.  Because  the  heat  stored  in  the  blades  can  be  carried  away  through  the 
cooling  air  convection,  the  blades  can  be  cooled  effectively,  as  shown  in  Fig.  1.1a.  But 
due  to  the  narrow  space  of  the  cooling  passage  at  the  blade  trailing  edge,  the  flow 
resistance  to  the  cooling  air  will  be  increased,  and  flow  rate  of  the  cooling  air  will  be 
reduced  there  accordingly.  As  a  result,  the  air  cooling  effect  at  the  trailing  edge  will  be 
decreased,  and  the  temperature  at  the  trailing  edge  will  be  substantially  higher  than  that 
of  other  regions  in  the  blade,  as  shown  in  Fig.  1.1b. 

(2)  Film  cooling,  in  which  the  cooling  air  not  only  passes  through  the  inner 
channels  in  the  blade,  but  also  flows  out  from  the  blade  surface  and  trailing  edge.  Air  film 
cooling  involves  the  injection  of  cooling  air  into  the  boundary  layer  of  high  temperature 
gas  on  the  blade  surfaces,  as  shown  in  Fig.  1.2.  Therefore,  it  provides  a  protective  layer  to 
the  blade  surfaces.  At  the  same  time,  the  temperature  gradient  in  the  cross-sectional  area 
of  the  blade  will  be  decreased  greatly.  But  the  air  film  cooling  not  only  needs  a  lot  of 
cooling  air,  but  also  adversely  affects  the  aerodynamic  and  thermodynamic  performance 
of  the  aerospace  turboengines.  For  example,  a  turbine  stage  with  the  film  cooling  would 
have  an  thermal  efficiency  of  about  85%.  But  for  a  turbine  stage  without  the  film  cooling, 
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the  thermal  efficiency  would  be  about  90%.  Therefore,  the  loss  in  thermal  efficiency  is 
substantial. 

(3)  Air  transpiration  cooling  or  porous  cooling,  which  involves  the  use  of  porous 
material.  The  cooling  air  in  this  case  is  forced  into  the  blade  surfaces  to  form  a  relatively 
cool  insulating  film,  as  shown  in  Fig.  1.3.  It  is  an  effective  cooling  method  for  the  blade, 
but  it  is  easy  to  lead  to  pore  blockage  due  to  material  oxidation.  It  would  also  cause  the 
blade  break  because  of  the  high  revolution  and  vibration  of  turboengines. 

The  third  approach  for  the  blade  cooling  is  the  heat  pipe  cooling  technology 
which  provides  a  potentially  effective  cooling  means  for  the  structure  exposed  to  high 
heat  fluxes.  Because  heat  pipes  have  operational  characteristics  of  having  a  very  high 
effective  thermal  conductance,  and  the  ability  to  act  as  a  thermal  flux  transformer  to 
create  an  isothermal  surface  of  low  thermal  impedance,  they  can  be  employed  in  the 
turbine  blades  as  a  cooling  device,  as  shown  in  Fig.  1.4. 
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cooling  air 


(a) 


(b) 


Fig.  1.2  Schematic  of  air  film  cooling  for  turboengine  blades.  (Bathie,  1984) 
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Fig.  1.3  Schematic  of  transpiration  air  cooling  for  turboengine  blades.  (Bathie,  1984) 
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1.2  Application  Backgrounds  and  Working  Principles  of  Radially  Rotating 
Miniature  Heat  Pipes 

Gas  turbine  blade  cooling  is  crucial  for  a  higher  combustion  temperature  and 
consequently  a  higher  engine  thermal  efficiency.  Air  cooling  technology  is  more 
convenient  and  effective,  but  the  coolant  flow  in  the  blades  affects  not  only  the  heat 
transfer  but  also  the  aerodynamic  and  thermodynamic  performance  of  the  aerospace 
turboengines.  At  the  same  time,  there  are  some  other  problems  associated  with  the  air 
cooling. 

1)  Air  cooling  effectiveness  is  reduced  towards  the  blade  tip  as  it  picks  up  heat. 
The  result  will  be  the  rise  of  the  general  level  of  blade  temperature  from  the  bottom  to  the 
tip.  Therefore,  it  will  result  in  a  large  temperature  gradient  between  the  blade  bottom  and 
the  tip. 

2)  Because  of  the  narrow  space  of  blade  trailing  edge,  the  flow  resistance  of  the 
cooling  air  will  be  increased,  and  flow  rate  of  cooling  air  will  be  reduced  there.  The 
temperature  gradient  over  the  cross-section  area  of  the  blade  will  be  increased,  as  shown 
in  Fig.  1.1b. 

3)  Air  cooling  technology  needs  about  2%  of  the  total  cooling  air  flow  from  the 
compressor  to  cool  the  high  temperature  components,  so  this  will  result  in  a  reduction  in 
the  efficiency  of  the  gas  turbine  equipment. 

In  order  to  reduce  the  temperature  gradient  in  turbine  blades  and  increase  the 
cooling  effectiveness,  combining  the  air  cooling  technology  with  radially  rotating 
miniature  heat  pipes  for  the  blade  cooling  is  an  innovative  approach.  As  mentioned 
earlier,  there  exists  a  large  temperature  gradient  over  the  cross-sectional  area  of  the  blade 
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and  between  the  tip  and  bottom  of  the  blade.  In  order  to  solve  this  problem,  a  radially 
rotating  miniature  heat  pipe  can  be  installed  in  the  blades,  as  shown  in  Fig.  1.5.  The 
radially  rotating  heat  pipe,  which  rotates  about  the  axis  of  revolution,  is  a  wickless  heat 
pipe  that  is  comprised  of  an  air-evacuated  hollow  container  and  an  amount  of  working 
fluid  sealed  within  the  container,  as  shown  in  the  Fig.  1.6.  The  centrifugal  force  generated 
by  the  rotating  motion  of  the  heat  pipe  is  utilized  for  the  return  of  the  liquid  condensate 
from  the  condenser  to  the  evaporator.  The  heat  in  the  blades  is  transferred  from  the 
evaporator  to  the  condenser,  where  the  heat  is  carried  away  by  the  cooling  air  convection. 
Because  the  radially  rotating  miniature  heat  pipe  has  a  very  high  thermal  conductance, 
the  temperature  gradients  over  the  cross-sectional  area  of  the  blade  and  between  the  tip 
and  bottom  of  the  blade  will  be  reduced  greatly  under  this  new  cooling  approach. 

In  addition,  the  top  portion  of  a  turbine  disk  will  operate  in  a  high  temperature 
condition  because  the  blades  of  a  gas  turbine  are  directly  attached  to  the  disk.  Therefore, 
the  radially  rotating  heat  pipe  can  be  also  applied  in  the  disc  cooling  application,  as 
shown  in  Fig.  1.7  (Cao,  1998). 
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TIP  CAP 


Fig.  1 .5  Schematic  of  air  convection  cooling  combined  with  radially  rotating  heat  pipe 
cooling  for  turboengine  blades. 


2.  ANALYSES  OF  RADIALLY  ROTATING  MINIATURE  HEAT 
PIPE  PERFORMANCE 

2.1  Analytical  Solution  for  the  Liquid  Film  Distribution  in  the  Condenser  Section 

As  mentioned  in  the  previous  chapter,  a  radially  rotating  heat  pipe,  which  rotates 
about  the  axis  of  revolution,  is  a  wickless  heat  pipe  that  is  comprised  of  an  air-evacuated 
hollow  container  and  an  amount  of  working  fluid  sealed  within  the  container.  Before  the 
heat  pipe  is  actually  applied  to  the  turbine  blade,  however,  the  fundamental  transport 
phenomena  related  to  various  heat  transfer  limitations  must  be  fully  understood. 

Fig.  2.1  schematically  illustrates  a  radially  rotating  heat  pipe  and  its  interior  working 
conditions.  The  assumptions  of  the  general  thin-film  condensation  analysis  are  as  follows: 

1 .  The  vapor  condensation  is  filmwise. 

2.  Liquid  subcooling  is  negligible. 

3.  Convective  effects  in  the  liquid  are  negligible. 

4.  The  liquid  film  thickness  is  much  smaller  than  the  vapor  space  radius. 

5.  Circumferential  velocity  and  temperature  gradients  are  negligible. 

6.  The  vapor  pressure  over  the  cross-sectional  area  of  the  miniature  heat  pipe  is 
uniform. 

7.  The  thermocapillary  flow  effects  are  negligible  in  the  presence  of  a  large  body 
force  field. 

However,  since  the  size  of  the  present  heat  pipe  is  very  small,  the  curvature  effect  of 
the  heat  pipe  wall  on  the  liquid  film  distribution  may  be  included.  Therefore,  a  cylindrical 
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coordinate  system  is  employed  in  the  present  study.  For  analyses,  a  cylindrical  differential 
control  volume  in  the  liquid  condensate  film  of  a  radially  rotating  heat  pipe  is 
schematically  shown  in  Fig.  2.2,  where  the  liquid  centrifugal  force  and  gravitational  force 

are  fc  =  2 7vp,co2  Zrdrdz  and  fg  =  2np, grdrdz  cos  6* ,  respectively.  A  static  force  balance  in 
z  direction  results  in  the  following  differential  equation: 


tyl.z 


dr 


p , .  2 7wdr  —  (pi .  h - —  dz)2xrdr  +  r, .  2nrdrdz  -  (r,  r  4 - —){2nr  +  dr)dz 

dz  ’*  dr 

+  Pi  (a2  Z  -  g  cos  Oyiizr dr  dz  =  0 


(2.1) 


where  Z  w  z0  +  z  is  the  revolving  radius  of  the  control  volume,  6  is  the  rotating  angle  of 
heat  pipe. 

Since  r»dr  and  2 tit  +  dr  «  2 nr ,  simplifying  Eq.  (2.1),  we  have 
dr,  „  dp, ,  . 

— —  +  — -- p,(co2Z  - gcos9)  =  0  (2.2) 

dr  dz 

The  static  force  balance  in  r  direction  is 

dp, , 

p,  r27trdz  —  (p,r  H - —dr)(27tr  +  dr)dz  =  0  (2.3) 

'  ’  dr 

Simplifying  the  equation  by  neglecting  higher  order  infinitesimals, 


dr 


=  0 


(2.4) 


Eq.  (2.4)  means  that  the  pressure  of  the  liquid  film  is  a  function  of  z  only.  The  shear  stress 
and  pressure  at  the  liquid  film  surface  can  be  studied  by  considering  the  momentum 
equations  in  both  z  and  r  directions  based  on  the  interface  control  volume  shown  in  Fig.2.3. 
We  get 
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—  dr, 

t,  R_s  2mdz  =  -Tr  R_s  2 nrdz  -  (w,  +  w,  R_s  )  —  2mdz  (2.5) 

dz 

—  dr, 

pr  nS  2  mdz  -  p,  R_g  (2m  +  dr)dz  +  (wr  +  w,  R_s  ) - 2  nrdz  -  0  (2.6) 

dz 

where  T,  is  the  liquid  mass  flow  rate  per  unit  circumferential  length,  wv  and  w,  R_s  are  the 
average  vapor  velocity  and  liquid  velocity  at  the  liquid  film  surface,  respectively,  r,  R_s 
and  r  R_s  are  the  shear  stresses  of  the  liquid  film  and  vapor  at  the  liquid-vapor  interface 

due  to  the  countercurrent  liquid  and  vapor  flows,  and  S  is  the  thickness  of  the  thin  liquid 
film,  as  shown  in  Fig.  2.3. 

dr. 

Because  -  is  very  small  and  can  be  ignored,  Eq.  (2.5)  and  (2.6)  can  be 

dz 

simplified  as  follows: 


Tl,R-S  ~  TvtR-S 

(2.7) 

Pl,RS  ”  P v,R-8 

(2.8) 

It  is  reasonable  to  assume  that  the  vapor  pressure  over  the  cross  section  of  the  heat  pipe  is 
uniform  since  the  cross-section  area  of  miniature  heat  pipes  is  very  small.  According  to 
Eqs.  (2.4)  and  (2.8), 

Pi(z)  =  PXz)  (2.9) 

This  means  that  the  pressures  of  vapor  and  liquid  in  the  r  direction  in  the  condenser 
section  are  the  same,  and  vary  only  in  z  direction.  Therefore,  Eq.  (2.2)  can  be  rewritten  as: 

^2Ll  +  ^Ei.-  pt(co2Z  -  gcosO)  =  0  (2.10) 

dr  dz 

The  boundary  conditions  for  the  above  equation  is 
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w,  =  0 


at  r  =  R 


(2.11) 


dw’i 

Ml  ~Z  =  ~Tl,R-g  =  Tr,RS 

or 


at  r  =  R  -  8 


(2.12) 


Integrating  Eq.  (2.10)  along  r  direction  yields  the  following  expression: 


-rlr  =  -^-r-  p,(co2Z-gcosO)r  +  c] 
dz 


(2.13) 


For  a  Newtonian  fluid, 


*7 ,  =  ~Mi 


(2.14) 


Plugging  Eq.  (2.14)  into  Eq.  (2.13)  and  integrating  it  along  r  direction  again  yield 


1  dp  y2C 

w  =  —  [— -  -  pt{co2Z  -  g  cos  0)]  —  +  —  r  +  c 
M,  dz  2  p, 


(2.15) 


According  to  boundary  conditions  (2.1 1)  and  (2.12),  we  have 
w,  =—[Ur2 - r2)-(R-r)(R-5)][p,(co2Z-gcos0 )-^rL] — -TVtR.s(R-r)  (2.16) 


Mi  2 


dz  p, 


Once  the  velocity  distribution  is  found,  the  liquid  mass  flow  rate  per  unit  circumferential 
length,  r, ,  is  evaluated  through  the  following  integral: 


m,  _  1  j* 

nD  2nR  •« 


^sp,w,2m-dr 


(2.17) 


For  the  radially  rotating  heat  pipe  with  a  high  speed,  the  liquid  film  thickness,  8 , 
is  much  smaller  than  the  heat  pipe  radius,  R.  Therefore,  r  and  R  are  roughly  equal  and  can 
be  canceled  in  Eq.  (2.17).  Substituting  Eq.  (2.16)  into  Eq.  (2.17)  and  integrating  it  yields 


T^^ip^Z-gcosO)- ~\ tvR_s 
p,  dz  3  pt  2 


(2.18) 
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The  following  relation  can  also  be  used  to  evaluate  the  mass  flow  rate  based  on  the  energy 
balance  by  neglecting  the  change  of  sensible  heat. 


r, 


(2.19) 


where  qc  and  qc  are  the  heat  flux  and  the  average  heat  flux  in  the  condenser  section, 
respectively.  Substituting  Eq.  (2.19)  into  Eq.  (2.18)  and  recasting  it,  the  solution  for  the 
liquid  film  thickness  is  found: 


~  P‘h-g.  [p, ( o2Z  —  g cos 0)  -  —]83  +  p,h*Ty-*  S2  + 1  =  0  (2.20) 

3/u,qcz  dz  2n,qcz 


Once  the  liquid  film  thickness  is  evaluated  from  Eq.  (20),  the  temperature  drop,  AT), 
across  the  liquid  film  can  be  calculated 


AT, 


R-8 

k, 


(2.21) 


The  local  heat  transfer  coefficient,  h, ,  and  local  Nusselt  number,  Nu:,  in  the 
condenser  section  are  defined,  respectively, 


h.  = 


<lc 


AT,  S 


,,  h.z  z 

Nu,  =  =  — 

*  k,  8 


where  k,  is  the  heat  conductivity  of  the  liquid  film. 


(2.22) 

(2.23) 
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2.2  General  Solution  for  the  Liquid  Film  Distribution  in  the  Condenser  Section  with 
a  Tilt  Angle 

The  basic  assumptions  of  the  mathematical  model  are  the  same  as  those  in  chapter 
2.1 .  In  general,  the  tilt  angle,  (f) ,  of  the  heat  pipe  should  be  close  to  90° ,  as  shown  in  Fig. 
2.4.  A  static  force  balance  in  the  z  direction  gives 

p,  2m/r  -  ( p ,  h - -  dz)2^rdr  +  z,  2 mdz  -  2n(r  +  dr)dz{z ,  +  — -  dr) 

dz  dr 

+  p,(o)2  Z  -  geos 8) sin </>2/trdrdz  =  0  (2.24) 

where  Z  =  z0  +  z  sin  <j>  is  the  revolving  radius  of  the  control  volume,  z0  is  the  revolving 
radius  of  the  coordinate  origin. 

Because  r»dr, 

2 n{r  +  dr)dz  =  2 m'dz  +  2 ndrdz  »  2 nrdz  (2.25) 

Plugging  Eq.  (2.25)  in  Eq.  (2.24)  yields 

p,  2 7trdr  -  (p,  +  dz)2wdr  +  z,  2nrdz  -  2nrdz(zl  H - -  dr) 

dz  dr 

+  p,{co2  Z  -  g cos 8) sin <j>27irdrdz  -  0  (2.26) 

Simplifying  Eq.  (2.26)  and  rearranging  it 

— -  =  — —  +  p  (co2Z  -  g  cos  8)  sin  (j>  (2.27) 

dr  dz 

A  static  force  balance  in  r  direction  would  give 

_ 

p,  2 7rrdz  -  ( pt  +  dr)2n(r  +  dr)dz  -  p,(co2  Z  -  g cos 8) cos <, t>2m~drdz  =  0  (2.28) 

dr 

Eq.  (2.28)  can  be  simplified  and  recast  as 
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=  -pt{co2Z  -  geos 9) cos (j> 
dr 


(2.29) 


Integrating  Eq.  (2.29)  from  r  to  R-S  yields 

Pi,r-s  ~ Pi  =  ~Pi(g>2 Z  - gcosO)cos0(r  -  R  +  S)  (2.30) 

Therefore  p,  =  pt  R_s  +  p,(o)2 Z  -  gcos9)cos0(r  -  R  + S)  (2.31) 

Differentiating  Eq.  (2.31)  with  respect  to  z  gives 

dPi_  _  Pi,r-6"  +  2  sjn^cos^(r  -p  +  £)  +  p^co2 Z  -  geos 0) cos <f>—  (2.32) 

dz  dz  dz 

Assuming  that  the  liquid  pressure  at  the  liquid-vapor  interface  is  equal  to  the  vapor 

pressure,  we  have 

P,*-s  =  Pr  (2.33) 

Substituting  Eq.  (2.33)  into  Eq.  (2.32)  gives 

^2.  =  ^Ez-+  piC0 2  sin <j) cos <j)(r  -  R  +  S)  +  p,(a>2 Z  -  geos 0) cos <f>—  (2.34) 

dz  dz  dz 

Plugging  Eq.  (2.34)  into  Eq.  (2.27),  we  have 

=  _£^x_  p)C02  sin  <j)  cos  <j>(r  -  R  +  5)  +  p,(co2  Z  -  geos 9)(sin0- cos 0—)  (2.35) 
dr  dz  dz 

.  dd  .  .. 

since  —  is  very  small, 
dz 


sm0-eos0 — «sin^  (2.36) 

dz 

Substituting  Eq.  (2.36)  into  Eq.  (2.35)  and  integrating  it  yield 

dv  i  r2 

t,  =  [ - -r  +  p,(o)  Z  -  geos 9) sin 0~\r  -  p,co  sin^cos^[ - (i?-<7)r]  +  c,  (2.37) 

dz  2 
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dw. 

Since  t,  =  —ju ,  — — 
dr 


replacing  r,  with  this  relation  in  Eq.  (2.37)  and  integrating  it. 


w,  =— !—  - p,(a>2 Z - g cos6)sin<f)]r2  +—a>2  sin^cos^[- - —  (7?-<5)/*2]-—  r  +  c, 

2  /u,  dz  /u,  6  2  /u,  ~ 

(2.38) 

The  two  constants  in  the  above  equation  can  be  determined  by  the  following  two 
boundary  conditions: 

w,  =  0  at  r=R  (2.39) 

_ 

ti,r-s  =  ~Tr,R-s  -  (w-  +  w,  R_s  )—L  at  r  =  R  -  8  (2.40) 


Since 


dr , 

dz 


is  very  small  and  negligible,  Eq.  (2.40)  can  be  simplified  as 


ti,r-s  =  -  tv,r-s  at  r  =  R  -  8  (2.41) 

After  determining  the  two  constants  in  Eq.  (2.38)  based  on  two  boundary  conditions,  the 
velocity  of  liquid  in  z  direction  is  obtained 

wl=—[^-pl(co2Z-gcos9)smmR-S){R-r)-UR1-r1)]-  —  TvR_s{R-r) 
jul  dz  2  //,  ' 


-  ^  (O1  sin  4, cos  </>[)- ( R  -  S)2 ( R  -  r )  -  )-  ( R  -  S)( R2  -  r 2 )  +  \  ( i?3  -  r 3 )]  (2.42) 


The  liquid  mass  flow  rate  per  unit  circumferential  length,  f, ,  is 


m  i  R 


R-6 


-^r  { [~- ~  P,  (a2  Z  ~  8 cos  0)  sin  fl(-  ±  RS3 
ju(R  dz  3 


23 


(2.43) 


-  p,(02  sin  (f> cos  RS*  -  — £5)} 

According  to  the  energy  balance,  f,  can  also  be  evaluated  by  the  following  relation 


Substituting  Eq.  (2.44)  into  Eq.  (2.43)  gives 


(2.44) 


hr  fitR  dz  3  24  2  j 


•l-*4  3  ,5, 


-  p,a>2  sin<z)cos^i?<54  £5)}  (2.45) 

Rearranging  Eq.  (2.45),  the  general  solution  for  the  liquid  film  thickness  is 


^3—  pfco2  sin^cos^5  + - /g^„  [p;(6)2Z-gcos^)sin^-^-]^4  -  fs^' 


40  p,Rqcz 


2Ap,Rqcz 


dz 


3p,Rqcz 


— y o,2a>2  sin^cos^<54 - fg^!  [p,(oj2 Z - g cos 0)  sin  ^  - ~L](^3  + — tv.rs^2  +1=0 

8  p,qcz  2p,qcz  dz  2  p,qcz 


(2.46) 


If  8  «  R,  plugging  r  s  R  into  Eq.  (2.43),  it  can  be  simplified  as 


r,=fD~=  I*** 


R-5 


c*  4  o  3  c*  2 

=  —  {-p,co2  sin  </>  cos  4>  —  -  [-^  -  P,  (CO 2  Z  -  g cos  0) sin  <f>\ —  -  r  R_s  — }  (2.47) 

pt  8  dz  3  2 

Substituting  Eq.(2.44)  into  Eq.  (2.47)  and  rearranging  it,  a  simplified  version  of  Eq. 

(2.46)  can  be  obtained 
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— — —  a)2  sin^cos^4 - [p,(a>2  Z  -  geos  6)  sin  <j>  -  3  + — —  rvR_sS2  +1=0 

8  Htqcz  3//;^cz  dz  2  n,qcz 

(2.48) 

Eq.  (2.48)  is  the  same  as  the  liquid  film  thickness  over  a  plane  geometry  with  a  tilt 
angle  (Cao  and  Chang,  1997),  since  the  cylindrical  surface  can  be  treated  as  a  plane 
surface  if  R»  8 .  Therefore,  Eq.  (2.48)  can  also  be  used  for  the  calculation  of  liquid  film 
thickness  for  a  plane  geometry  with  a  tilt  angle. 

If  the  tilt  angle,  <f> ,  is  90° ,  then  the  equation  for  the  liquid  film  thickness  can  be 
simplified  further  as 

- fg^!  [p^ai1  Z  -  geos  0)  sin  -^-]8 2  ■+ — tv r_s52  +1  =  0  (2.49) 

2>^qcz  dz  2  n,qez 

Eq.  (2.49)  and  Eq.(2.20)  are  exactly  the  same  because  both  of  them  are  the 
equation  of  the  liquid  film  thickness  for  a  cylindrical  heat  pipe  with  <j>  =  90° . 

Comparing  Eq.  (2.48)  with  Eq.  (2.46),  it  is  clear  that  Eq.  (2.48)  can  be  obtained 
from  Eq.  (2.46)  if  the  first,  second  and  third  terms  on  the  left-hand  side  of  Eq.  (2.46)  are 

dropped.  Also,  Eq.  (2.49)  is  obtained  from  Eq.  (2.48)  with  </>  =  90° .  For  these  reasons, 
Eqs.  (2.48),  (2.20)  and  (2.49)  are  three  special  cases  of  Eq.  (2.46).  In  other  words,  Eq. 
(2.46)  is  the  general  relation  for  liquid  film  distributions  in  radially  rotating  heat  pipes. 

After  the  liquid  film  thickness,  8,  is  obtained  from  the  equations  for  the  liquid 
film  thickness,  the  temperature  drop  across  liquid  film  thickness  is 

Rqc  In — — — 

AT)  =Tg-Tp= - (2.50) 

/C/ 
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If  the  temperature  distribution  of  liquid  film  is  linear,  the  temperature  drop  across  the 
liquid  film  can  be  expressed  as 


A  T 


The  local  heat  transfer  coefficient,  h, ,  and  local  Nusselt  number,  Nu, , 
section  can  be  the  same  as  Eqs.  (2.22)  and  (2.23),  respectively. 


(2.51) 

in  the  condenser 
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Fig.  2.4  Schematic  of  a  differential  control  volume  in  the  condensate  film. 


2.3  Axial  Vapor  Flow  Analyses  along  the  Heat  Pipe  Length 


In  the  foregoing  analytical  solutions  for  the  liquid  film  distributions  on  the  heat 
pipe  wall,  both  the  shear  stress  at  the  liquid-vapor  interface  and  the  vapor  pressure  gradient 
along  the  heat  pipe  length  are  present.  To  complete  the  analysis,  solutions  for  the  vapor 
flow  must  be  obtained.  Also,  as  mentioned  in  the  introduction,  the  centrifugal  force  may 
play  an  important  role  in  the  vapor  temperature  drop  along  the  heat  pipe  length.  In  some 
extreme  cases,  this  temperature  drop  may  present  a  heat  transfer  limitation  for  the  heat  pipe 
performance.  As  a  result,  it  is  necessary  to  obtain  the  temperature  drop  along  the  heat  pipe 
length.  Consider  a  one-dimensional,  laminar,  and  incompressible  vapor  flow  along  the  heat 
pipe  length  under  the  steady-state  operating  condition,  as  shown  in  Fig.  2.5,  where 

fc  =  p  a,2  ZnR]  sin^iz  is  the  centrifugal  force  of  the  vapor,  and 
/  =  pxgnR]  sin  <j> cos  Gdz  is  the  gravitational  force  of  the  vapor.  The  momentum  equation 
for  the  vapor  flow  is 

pvnR)  sin^  -  ( pv  +  ^Ldz)nR]  sin^  +  rv  R  2 xRvdz  +  pv(co2  Z  -gcosG)xR2  sin <jxiz 
dz 


=  mv  wv  —  [mv  w,.  +  — (mv  w^dz\ 
dz 


Eq.  (2.52)  can  be  simplified  as 


dPr  2  ^ 

dz  Rv  sin  <j)  ’’’ 


^  ^  2  v  „ — „  n\  d  ...  \  1 


tv,r„  +  Pi  (®  z  -  S  cos  G)  -  — ■  (m  v  Wv ) —J 


nR2  sin (j> 


Since  the  shear  stress  of  the  vapor  flow  is 


(2.52) 


(2.53) 
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(2.54) 


1  -2  2  wrRrp  //,,  1  -2  jUvWv 

^  =-fprWy  =/( - —  Xr=-r — =(/Re)r  — 

2  //,.  2wrRrpv  2  4i?v 


where  /?,,  =  R  -  8  is  the  vapor  space  radius  in  the  heat  pipe,  mr  is  the  vapor  mass  flow 

rate,  wv  is  the  average  vapor  velocity  over  the  heat  pipe  cross-sectional  area,  and  /  is  the 
laminar  vapor  skin-friction  coefficient. 

A  thermal  balance  in  the  condenser  section,  as  shown  in  Fig.  2.6,  gives 


2 nRvqcz  2  - 

mv  = - —  =  7rRr  px  w,. 

h* 


(2.55) 


Therefore,  w,.  = 


2  g^_ 
Rvpvh 


.  Substituting  it  into  Eq.  (2.54)  yields 


Is 


Since 


tv.r..  =(/Re)v 


Mv<lcz 

2  Rlpvhfg 


d(mv  w>.) 
dz 


_  -"2 
2  _  —  dwy  8 nqc  z 


=  2nRv  pv  w,. 


dz  Pyh/g 


Plugging  Eqs.  (2.56)  and  (2.57)  into  Eq.  (2.53)  gives 


(2.56) 


(2.57) 


dPv  =  (/Re),.  ^'?£~  +  Pi(a>2 Z - gcosQ) -  8<7c  2 


dz  . Rlpyhfgsm<j>  ~  '  pvRyhjg  sin^ 


(2.58) 


In  general,  the  above  equations  are  valid  when  the  following  two  conditions  in  terms  of  the 
vapor  Reynolds  number  and  Mach  number  are  satisfied: 

2  RyQ 


Re.,  = 


’’  Avjuvh 


<2300 


(2.59) 


7s 
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Mr  = - <  03 

AyPyhfg^y  rRgTy  (2.60) 

where  Ar  is  the  cross-sectional  area  of  the  vapor  space,  yr  is  the  vapor  specific  heat  ratio, 
Rg  is  the  gas  constant  of  the  vapor  (Chi,  1976).  These  two  conditions  are  usually  satisfied 
for  high-temperature  heat  pipes  having  a  small  diameter  and  operating  at  a  steady  state. 

Assuming  that  the  heat  addition  along  the  evaporator  and  the  heat  removal  along 
the  condenser  are  constant,  the  axial  mass  flow  rate  of  vapor  can  be  found  by  a  mass  and 
energy  balance  under  the  steady  state  (Faghri,  1995). 


where  q,  is  axial  heat  transfer  rate  of  the  vapor. 

Substituting  Eq.  (2.61)  into  Eq.  (2.55)  and  rearranging  it 

0  <z<Lc 

Lc<z<Lc+La  (2.62) 

Lc  +  L  a  —  z  —  L 

Plugging  Eq.  (2.62)  into  Eq.  (2.58)  and  integrating  it  from  0  to  L,  the  vapor  pressure  drop 
along  the  heat  pipe  is 


r 

QcZ 

2tcRxLc 

m"  hfg 

Qc  K 

2  t:Rv 

Lch/g  2ttRv 

Qc 

hfg  L  —  z 

h/g 

2nRx  L-(LC+La)hfg 
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AP’  =  f 


dp* 

dz 


dz 


f  (/  Re>.  #^77-  &  +  f  P,  Z  -  g  oos  0)dz  -  -±-  f  -%£-  * 

-1  R*P*K  *  sin^  •*>  pvR;h~fg 


sin  (j) 


(2.63) 

Since  the  vapor  momentum  are  very  small,  and  can  be  negligible,  the  vapor  pressure  drop 
along  heat  pipe  can  be  simplified  as 


1  _  ,Le  +2La  +4)  +  ^ vO)2l(z0  +  -|sin^)-p rLgcos9 


Apr=—FVQC( 

sin^  2 


1  2  - 
-r-7  FyLeJfQc  +  A®  LZ°  -pyLgcosO 
sin^z) 


(2.64) 


where  L  is  the  heat  pipe  length,  LeJf  =  ^(Ic  +2La  +  Le)  is  the  effective  length  of  the 


ff  jj  _  J 

heat  pipe,  Fv  = - —  -  -  is  vapor  frictional  coefficient,  Za  =  z0  +  —  L  smcj)  is  average 


2  nRvpvh 


fg 


revolving  radius  of  the  heat  pipe,  and  Qc  =  Qe  =  nDLcqc  is  the  total  heat  transfer  rate  at 
the  condenser  or  evaporator  section  at  a  steady  state. 

As  mentioned  earlier,  the  vapor  Reynolds  number  is  in  general  less  than  2300.  But 
with  an  increase  in  the  total  heat  input  Q,  the  Mach  number  may  be  greater  than  0.3  when 
the  heat  pipe  operating  temperature  is  relatively  low.  In  this  situation,  the  vapor 
compressibility  should  be  taken  into  consideration.  Based  upon  the  bulk  fluid  properties, 
the  ratio  of  the  drag  coefficient  for  a  compressible  flow,  frc ,  to  that  for  an  incompressible 


flow,  /,, . ,  at  the  same  vapor  Reynolds  number  can  be  corrected  by  the  following  equation 
(Chi,  1976): 
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(2.65) 


f \ v,c  Y V  1 

—  =  (1  +  — — Mv)  2 

Jr,i 

Substituting  Eq.  (2.65)  into  the  relation  for  the  vapor  frictional  coefficient,  a  modified 
vapor  friction  coefficient  is  obtained: 


F„  = 


(/Re,. K  „  ,r,.+l1/2 


(1  + 


m;T 


(2.66) 


2nRlpyhfs  v‘  2 

In  fact,  the  vapor  pressure  drop  caused  by  the  gravity  is  very  small,  and  can  be 
negligible.  Therefore,  the  vapor  pressure  drop  along  the  heat  pipe  length  becomes 


A p,=^LQ,*PyLZ. 

sin  (p 


(2.67) 


The  Clapeyron  equation  that  relates  the  saturated  temperature  to  the  saturated  pressure  is 

hJ*P* 


A Pv  = 


T 


■AT 


(2.68) 


Substituting  Eq.  (2.68)  into  Eq.  (2.67),  the  vapor  temperature  drop  along  the  heat  pipe 
length  is 

T 


TyL  ffF  T  ,  - 

AT-  =  Ap,  =  « ;--Q,+-pa'LZ.  =  AT  ,+ AT 

hfgP v  pyh/gsm<p  hfg 

The  dimensionless  temperature  drop  of  the  vapor  along  the  heat  pipe  length  is 


(2.69) 


A  Tv  LeffFv 


Tv  Pvhfgsm<p  h 


k  =  Ar,  +  Ar. 


(2.70) 


Js 


where  A Tv  is  the  total  vapor  temperature  drop  along  the  heat  pipe  length,  Tr  is  the 
operating  temperature  of  the  heat  pipe,  A Tf  is  the  dimensionless  temperature  drop  due  to 
friction  at  the  liquid-vapor  interface,  and  A Tc  is  the  dimensionless  temperature  drop  due 
to  the  vapor  centrifugal  force. 
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Fig.  2.6  Schematic  of  an  overall  thermal  balance  in  the  condenser  section. 


2.4  Analytical  Results  for  the  Temperature  Drops  and  Heat  Transfer  Limitations 


Temperature  drops  across  the  liquid  film  or  along  the  heat  pipe  length  are  very 
important  criteria  for  designing  radially  rotating  miniature  heat  pipes  with  a  high 
operating  temperature  and  rotating  speed.  In  order  to  prove  the  feasibility  of  the  turbine 
blade  cooling  application,  analytical  results  are  presented  for  radially  rotating  miniature 
heat  pipes  working  under  high  heat  flux  and  high  rotating  speed  conditions.  The  ranges  of 
heat  inputs,  geometric  dimensions,  and  the  rotating  speeds  are  as  follows:  40  W  <  Q  < 

200  W,  1  mm  <di  <  8  mm,  6.26  <co2Za  / g  <  6.26  xlO5,  L  =  80  mm,  Lc  =  Le=  40  mm, 

La  -  0 ,  and  (j)  =  90° .  It  should  be  pointed  out  that  these  ranges  are  only  close  to  the 
turbine  blade  cooling  conditions,  not  real  design  values  (Harman,  1981).  It  is  clear  that, 
according  to  Eqs.  (2.46)  through  (2.70),  the  liquid  film  thickness,  temperature  drop  across 
the  liquid  film,  and  vapor  temperature  drop  along  the  heat  pipe  length  are  mainly  the 
function  of  shear  stresses,  heat  fluxes,  liquid  and  vapor  thermophysical  properties,  liquid 
and  vapor  centrifugal  forces,  heat  pipe  sizes,  and  heat  pipe  operating  temperatures.  The 
operation  of  a  heat  pipe  depends  on  the  heat  load  in  the  evaporator  and  the  cooling 
condition  in  the  condenser.  In  general,  the  operating  temperature,  Tv ,  of  a  heat  pipe  can 
be  controlled  by  the  flow  rate  of  cooling  air  in  the  condenser  section  and  heat  input  in  the 
evaporator  section.  The  operating  temperature,  Tr ,  will  be  going  up  if  the  flow  rate  of  the 
cooling  air  is  reduced  or  the  heat  input  is  increased. 

Fig.  2.7  illustrates  the  maximum  temperature  drops  across  the  liquid  film  versus 

the  dimensionless  centrifugal  forces,  co2Zal  g,  for  sodium  as  the  working  fluid.  The 
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working  temperature,  Ty ,  and  the  total  heat  input,  Q,  are  kept  as  constants  at  1100K  and 
200  W,  respectively.  The  dimensionless  geometric  parameter,  d,l  L,  is  treated  as  a 
variable.  For  the  results  presented  in  Fig.  2.8,  however,  dt  I L  is  kept  as  a  constant,  and 
the  operating  temperature  is  considered  to  be  a  variable.  It  is  observed  from  Figs.  2.7  and 

2.8  that  at  a  relatively  high  rotating  speed  or  large  revolving  radius,  co2  Z„  / g>  150,  the 
maximum  temperature  drop  across  the  liquid  film  is  negligible.  However,  a  small  heat 
pipe  diameter  or  a  low  heat  pipe  operating  temperature  may  increase  the  temperature 
drop.  Still,  under  normal  turbine  working  conditions,  the  maximum  temperature  drop  is 
much  less  than  1  K.  This  means  that  the  commonly  encountered  condenser  limitation  for 
low-temperature  heat  pipes  would  never  be  a  serious  problem  for  the  high-temperature 
heat  pipes  studied  in  this  report. 

Calculations  are  then  made  for  the  vapor  temperature  drops  along  the  heat  pipe 
length  with  sodium  and  potassium  as  the  working  fluids.  Fig.  2.9  shows  vapor 
temperature  drops  as  a  function  of  dimensionless  centrifugal  forces  at  different  operating 
temperatures.  The  dimensionless  heat  pipe  geometric  parameter  and  the  total  heat  input 
are  kept  as  constants  (d,l L=  0.025  and  Q  =  200  W).  It  is  clear  from  the  figure  that  the 
vapor  temperature  drops  increase  with  an  increase  in  the  centrifugal  force  or  in  the 
operating  temperature.  At  the  same  time,  it  can  be  seen  that  a  heat  pipe  with  sodium  as 
the  working  fluid  would  have  a  smaller  vapor  temperature  drop  than  the  heat  pipe  with 
potassium  as  the  working  fluid.  For  many  cases,  the  vapor  temperature  drop  is  on  the 
order  of  10-40  K.  This  temperature  drop  is  considerably  high  for  a  low-temperature  heat 
pipe.  However,  for  the  high-temperature  heat  pipes  studied  in  this  report,  this  temperature 
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drop  is  relatively  small  compared  to  the  average  heat  pipe  operating  temperature  level. 

The  vapor  temperature  drops  as  a  function  of  heat  pipe  operating  temperatures 
and  dimensionless  size  are  shown  in  Fig.  2.10.  At  a  given  heat  pipe  operating 
temperature,  the  vapor  temperature  drop  is  increased  when  the  heat  pipe  diameter  is 
decreased  due  to  a  higher  friction  at  the  liquid-vapor  interface.  If  the  heat  pipe  operating 
temperature  is  less  than  1100K,  the  vapor  temperature  drops  for  the  miniature  heat  pipe 
with  dj  /  L  =0.0125  would  increase  sharply.  However,  for  a  heat  pipe  with  a  dj  L 

greater  than  0.02  or  d:  greater  than  1.5  mm,  the  influence  of  the  operating  temperature 
level  on  the  vapor  temperature  drop  is  relatively  small. 

Fig.  2.11  illustrates  the  influence  of  the  heat  input  on  the  vapor  temperature  drop 
along  the  heat  pipe  length.  In  general,  the  influence  of  the  heat  input  on  the  vapor 
temperature  drop  is  relatively  small  when  the  operating  temperature  level  of  the  heat  pipe 
is  relatively  high.  However,  for  a  heat  pipe  operating  at  a  relatively  low  temperature 
level,  the  vapor  temperature  drop  is  substantially  increased  when  the  heat  input  is  high. 

According  to  Eq.  (2.70),  the  total  vapor  temperature  drop  is  due  to  the 
contributions  of  friction  at  the  liquid-vapor  interface  and  the  vapor  centrifugal  force  in 
the  heat  pipe.  The  contribution  of  the  centrifugal  force  is  illustrated  in  Fig.  2.12  as  the 
ratio  of  the  vapor  temperature  drop  due  to  the  vapor  centrifugal  force  to  the  total  vapor 
temperature  drop.  At  a  low  rotating  speed,  the  contribution  of  the  centrifugal  force  is 
relatively  small,  but  at  a  high  rotating  speed,  the  contribution  of  the  centrifugal  force 
becomes  a  dominant  factor.  This  is  especially  obvious  for  a  heat  pipe  having  a  relatively 
large  diameter. 
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As  illustrated  by  the  foregoing  analytical  results,  the  vapor  temperature  drop 
along  the  heat  pipe  length  may  present  a  heat  transfer  limitation  when  the  heat  pipe  is  of 
very  small  diameter,  with  a  relatively  low  operating  temperature  level  and  high  rotating 
speed.  However,  for  a  miniature  heat  pipe  having  a  diameter  on  the  order  of  1 .5  to  2  mm, 
the  heat  transfer  limitation  would  not  be  encountered  at  normal  rotating  speeds  and 
operating  temperature  levels.  This  indicates  that  the  radially  rotating  miniature  heat  pipe 
is  feasible  for  the  turbine  blade  cooling  application. 

The  entrainment  limit  is  due  to  the  interaction  between  the  countercurrent  liquid 
and  vapor  flows  occurring  at  the  liquid- vapor  interface.  This  limit  may  be  a  major 
concern  when  the  heat  input  is  high  and  the  heat  pipe  dimensionless  sizes  are  small.  In 
this  case,  the  shear  stresses  at  the  liquid-vapor  interface  retard  the  return  of  the  liquid 
from  the  condenser  to  the  evaporator  section.  The  heat  pipe  is  said  to  have  reached  the 
entrainment  limit  when  the  liquid  is  prohibited  from  returning  to  the  evaporator  section. 
A  comprehensive  review  regarding  the  entrainment  limit  for  two-phase  closed 
thermosyphons  is  given  by  Faghri  (1995)  and  Peterson  (1994).  Various  correlations  for 
the  evaluation  of  this  limit  are  also  included  in  the  literature.  For  current  wickless  rotating 
heat  pipes,  the  mechanism  of  the  entrainment  limit  should  be  similar  to  that  of 
thermosyphons  with  the  gravitational  force  being  replaced  by  the  rotating  centrifugal 
force.  Fig.  2.13  shows  the  entrainment  limits  for  the  same  heat  pipe  studied  earlier  at 
different  rotating  speeds.  The  calculation  was  made  based  on  Wallis’  correlation  (Faghri, 
1995) 


Qeiu 
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and  the  correlation  from  Faghri  et  al.  (1989) 

%■  =  (p,/p,)° "  (tanh!  Boi )h/g[(r0  +  z)oV(A  -  p,)tW‘  +  p?‘ ]'!  (2.72) 

for  sodium  and  potassium  as  the  working  fluids.  For  the  rotating  heat  pipe  with  a  small 
dimensionless  size,  dt  /  L=  0.0125,  the  entrainment  limit  is  very  sensitive  to  the  rotating 
speed,  and  it  may  be  a  concern  when  the  rotating  speed  is  very  low.  However,  when  the 
rotating  speed  is  moderately  high,  the  entrainment  limit  is  increased  to  more  than  100 
W/cm2,  which  is  acceptable  for  turbomachinery  applications.  The  values  of  the 
entrainment  limit  that  are  calculated  by  the  Wallis’  correlation  are  close  to  those  from 
Faghri  et  al.  (1989),  when  the  rotating  speed  is  low.  However,  at  high  rotating  speeds, 
large  deviation  occurs  with  both  working  fluids.  A  similar  trend  was  observed  also  by 
Faghri  (1995)  for  thermosyphons. 
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Fig.  2.7  Maximum  temperature  drops  across  the  liquid  film  for 
sodium  as  the  working  fluid  (T  =  1 100  K,  Q  =200  W). 
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Fig.  2.8  Maximum  temperature  drops  across  the  liquid  film  for  sodium 
as  the  working  fluid  ( djL  =  0.025 ,  Q  =200  W). 
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(b)  Potassium 


Fig.  2.9  Vapor  temperature  drops  along  the  heat  pipe  length  at  different 
rotating  speeds  {dJ L  =0.025,  Q= 200  W). 


(b)  Potassium 


Fig.  2.10  Vapor  temperature  drops  along  the  heat  pipe  length  at  different 
heat  pipe  dimensionless  sizes  (Q= 200  W,  co2  Zaj g  =  6.26  x  104 ). 
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Vapor  temperature  drops  along  the  heat  pipe  length  at  different 
(dJL  =0.025,  (o2Za!g  =  3.91xl03). 


C02Za  /  g 

(a)  Sodium 


0)2Za  /  g 
(b)  Potassium 

Fig.  2.12  Ratio  of  the  vapor  temperature  drops  due  to  the  centrifugal  force  to  the 
total  temperature  drop  ( Ty =1 1 00  K,  £>=200  W). 
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Fig.  2.13  Entrainment  limitations  as  a  function  of  rotating  speeds 
(Tref- 937  K  for  potassium  and  TKf  =1073  K  for  sodium). 
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3.  TEMPERATURE  DISTRIBUTIONS  ALONG  THE  HEAT  PIPE 
AND  DIFFUSE  EFFECTS  OF  NON-CONDENSABLE  GASES 


3.1  Analysis  of  Vapor  Temperature  Distribution  along  the  Heat  Pipe 


Temperature  distributions  along  the  heat  pipe  length  play  an  important  role  in  the  heat 
transfer  efficiency  of  the  heat  pipe.  Vapor  temperature  drop  for  traditional  heat  pipes 
without  rotation  are  mainly  due  to  the  vapor  friction.  The  vapor  temperature  is  almost 
uniform  since  the  vapor  friction  is  generally  small.  But  for  radially  rotating  heat  pipes, 
the  temperature  drop  is  caused  not  only  by  the  vapor  friction  but  also  by  the  centrifugal 
force.  Therefore,  the  temperature  drop  for  a  radially  rotating  heat  pipe  is  larger  than  that 
for  a  traditional  heat  pipe. 

For  the  vapor  momentum  equation  in  chapter  2.3,  the  last  term  of  the  right-side  in  Eq. 
(2.53)  is  usually  negligible  compared  to  the  other  terms,  and  the  pressure  gradient  Eq. 
(2.53)  becomes 


dpv 

dz 


=  (/  Re)„ 


cz 


Rvprhfgsm<j) 


+  p^ccdZ  -  gcosd) 


(3.1) 


In  addition,  for  the  practical  application  of  blade  cooling,  the  adiabatic  length,  La , 
is  zero.  Therefore,  the  heat  pipe  is  divided  only  into  the  evaporator  and  condenser 
sections.  Integrating  Eq.  (3.1)  from  0  to  z  based  on -the  mass  and  energy  balance  under 
the  steady  state  in  Eq.  (  2.62)  and  La  =  0  yields, 
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sirup  '  2 nR*pvhfg  *  Lc 
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2sirup  L-Lc  2 

The  vapor  gravitational  force  can  be  neglected  compared  to  the  centrifugal  force  in 
Eqs.(3.2)  and  (3.3).  According  to  Clapeyron  Equation,  we  have 


A  r  = 
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since 
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substituting  Eqs.  (3.4)  and  (3.6)  into  (3.5)  and  rearranging  it  yield 


Tr.  =  7/„  +  ATr.  = 
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7? 


(3.7) 


where  ATr  is  the  total  vapor  temperature  drop  along  the  heat  pipe  length,  7/o  and  Tv  I 

are  the  temperatures  at  the  condenser  and  evaporator  ends,  respectively. 

Eq.  (3.7)  gives  an  analytical  solution  for  the  temperature  distributions  along  the 
heat  pipe  length.  It  should  be  pointed  out  that  Eq.  (3.7)  only  considers  the  effects  of  vapor 
friction  and  centrifugal  force  on  the  temperature  distribution.  If  there  are  some  non¬ 
condensable  gases  in  the  heat  pipe,  the  diffuse  effects  of  non-condensable  gases  on  the 
temperature  distribution  must  be  included. 


49 


3.2  Analytical  Solution  of  Temperature  Distribution  along  the  Heat  Pipe  Length 

Including  the  Diffuse  Effects  of  Non-Condensable  Gases 

Because  of  the  wickless  structure  and  a  very  small  volume  of  the  radially  rotating 
miniature  heat  pipes,  it  is  possible  to  have  a  small  amount  of  non-condensable  gases 
remaining  in  the  heat  pipe  during  the  procedure  of  charging  and  sealing.  Non¬ 
condensable  gases  will  be  driven  to  the  condenser  end  by  the  vapor  from  the  evaporator, 
which  prevent  the  vapor  to  fill  the  entire  condenser  section.  Therefore,  the  temperature 
distribution  along  the  condenser  wall  will  not  be  uniform  any  more,  and  a  large 
temperature  drop  will  appear  in  the  condenser  section  as  a  result  of  the  existence  of  non¬ 
condensable  gases.  Marcus  (1971)  first  developed  a  numerical  model  based  on  one¬ 
dimensional  diffusion  of  the  vapor-gas  system  coupled  with  one-dimensional  axial  heat 
conduction  in  the  wall  for  variable  conductance  heat  pipes.  Numerical  results  indicated 
that  the  axial  heat  conduction  plays  a  dominant  role  in  the  determination  of  the 
temperature  distribution  along  the  heat  pipe  wall. 

For  a  radially  rotating  heat  pipe,  the  vapor-gas  interface  between  the  vapor  and  non¬ 
condensable  gases  is  assumed  to  be  a  flat  front.  Fig.  3.1  shows  schematically  a  thermal 
balance  for  the  heat  pipe  wall  in  the  condenser  section.  A  control  volume  of  the  cross- 

sectional  area,  ;r(r02  -  r,2 ) ,  and  infinitesimal  length,  dz,  is  considered.  Based  on  the 
conservation  of  energy, 

n(ro  ~ r? )kp  -j-  +  2xr0hc (T  -T,)dz  =  ntf  -  r?  )k  dz)  +  £)  dz  (3.8) 

dz  dz  dz 

Eq.  (3.8)  can  be  simplified  as 
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(3.9) 


-  r,  )*,  4  -  2*r<A (T„-T,)  +  Q' =  0 

uZ 

Since  the  radial  temperature  drop  in  the  heat  pipe  wall  is  much  less  than  the  axial 
temperature  drop,  the  radial  temperature  drop  can  be  neglected,  that  is,  the  temperature  in 
the  radial  direction  of  heat  pipe  wall  is  uniform,  Tp=T .  Eq.  (3.9)  can  be  rearranged  as 


d2T 
dz 2 


-m2(T-Ts- 


Q 

2  7tr0hc 


)  =  0 


(3.10) 


where 


m2  = 


2r0hc 


kp(r0  -r,  ) 


(3.11) 


hc  is  the  mean  heat  transfer  coefficient,  Q  is  the  heat  output  in  the  condenser  section 
without  non-condensable  gases,  m  is  a  coefficient,  T  is  the  axial  temperature  which  is  a 
function  of  2,  Tp  and  Ts  are  the  outer  surface  temperature  of  the  heat  pipe  and 
surrounding  temperature,  respectively. 

Assuming  the  distribution  of  Q  is  uniform  in  the  portion  of  the  condenser  section 
without  non-condensable  gases,  and  zero  in  the  portion  of  the  condenser  section  with 
non-condensable  gases, 


o 

II 

O) 

for  0 <z<Lcn 

(3.12) 

Q  = 

4-4, 

for  Lcn<z<Lc 

(3.13) 

where  Qc  is  total  heat  output  in  the  condenser  section;  Lc  is  the  lengths  of  the  total 
condenser  section  and  Lc  n  is  the  portion  with  the  non-condensable  gases,  respectively. 

At  the  inlet  of  the  condenser  section,  the  axial  heat  conduction  is  negligible, 


51 


a=2Sr,(4-i„)J,(rIi- 

~TS) 

(3.14) 

where  TLc  is  the  temperature  at  the  condenser  inlet. 

Substituting  Eq.  (3.14)  into  (3.13)  and  recasting  it 

Q  =T  j 

l7rr0hc 

at  L  <  z  <  L 

(3.15) 

plugging  Eqs.  (3.15)  and  (3.12)  into  (3.10)  yields 

d  T  2  /rp  rp  k  n 

m  (T  T,)  =  0 
dz 

at  0  <  z  <  L 

c,n 

(3.16) 

d  T  2  ^  n 

at  L  <  z  <  L 

c,n  c 

(3.17) 

For  a  radially  rotating  miniature  heat  pipe,  the  boundary  conditions  are 

dz 

at  z  =  -oo 

(3.18) 

ii 

at  z  =  Lc 

(3.19) 

According  to  the  balance  of  total  thermal  energy  in  the  condenser  section. 

l’2m-MT-Tc)ctz  =  Q, 

(3.20) 

Solving  Eqs.  (3.16)  and  (3.17)  yields 

T  =  a  exp (mz)  +  b  exp (~mz)  +  Ts 

at  0  <  z  <  Lc  n 

(3.21) 

T  -  c  exp  (mz)  +  d  exp  (-mz)  +  TLc 

at  Lc  n  <  z  <  Lc 

(3.22) 

Applying  boundary  conditions  Eqs.  (3.18)  and  (3.19)  to  determine  the  constants  in  Eqs. 
(3.21)  and  (3.22)  gives 

b  =  0,  and  d  =  -c  exp(2 mLc ) 
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Substituting  coefficients,  b  and  d,  into  Eqs.  (3.21)  and  (3.22)  results  in 

T  =  aexp(mz)  +  Ts  at  0  <  z  <  Lcn  (3.23) 

T  =  c{cxp(mz) - exp[/w(2Zc  - z)]}  +  at  Lcn  <z<Lc  (3.24) 

Because  of  the  continuity  of  axial  temperature  distribution  along  the  heat  pipe  length, 
the  temperature  in  Eqs.  (3.23)  and  (3.24)  must  be  equal  at  z  =  Lc  n , 


a  exp (mLc  n )  +  T,=  c{txp(mLcn )  -  exp[m(2Lc  -  Lc  n )] }  +  TK 


(3.25) 


Therefore, 


c  = 


aexp(mLcn)-(T,c  -Ts) 
exp (mLc  n )  -  exp[m(2Lc  -  Lc  n )] 


(3.26) 


Substituting  Eqs.  (3.23)  and  (3.24)  into  (3.20)  and  rearranging  it  give 


a  f'  exp(mz)£fe  +  f'  {(7^  -Ts)  +  c[exp(mz)  - exp(m(2Zc  -  z))]}dz  =  - . —  -  (3.27) 

•o  * c.11  c  2nr0hc 

Plugging  Eqs.  (3.14)  and  (3.26)  into  Eq.  (3.27),  integrating  and  recasting  it  yield 

(Tl  -  Tx){2 exp(mLc )  -  exp {mL  )-  exp[m(2Zc  -  Zc  „)]} 

a  =  — -  (J.zo) 

2  exp[m(Zc  +  Lc  n )]  +  Qxp[m(2Lc  -  Lc  n )]  -  2  exp(2wlc )  -  exp (mLc  n ) 

Substituting  the  coefficients,  a  and  c,  into  Eqs.  (3.23)  and  (3.24)  gives  the  temperature 

distribution  along  the  condenser  section  under  the  diffuse  effects  of  non-condensable 


gases. 


If  Lcn  =  0 ,  there  is  no  non-condensable  gas  in  the  condenser  section.  Substituting 


Lcn  =  0  into  Eqs.  (3.28)  and  (3.26),  we  get 

a  =  T,c  -  Ts  and  c=0 

Plugging  a  =  Tl  —Ts  and  c=0  in  Eqs.  (3.23)  and  (3.24)  yields 
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0  <z<L 


T  =  T, 

The  above  relation  means  that  the  temperature  distribution  along  the  condenser 
section  is  uniform  under  the  condition  of  no  temperature  drops  caused  by  the  friction  and 
centrifugal  force.  At  the  same  time,  it  also  proves  the  conclusion  in  chapter  2,  in  which 
the  temperature  drops  along  the  heat  pipe  are  mainly  caused  by  the  friction  and 
centrifugal  force  if  there  are  no  non-condensable  gases  existing  in  the  heat  pipe. 

It  should  be  pointed  out  that  since  the  radial  temperature  drop  in  the  heat  pipe 
wall  is  much  smaller  than  the  axial  one  and  is  neglected  in  the  governing  equation  (3.10), 
the  temperature  in  the  heat  pipe  wall  can  represent  the  vapor  temperature  in  this  region 
proximately. 

Substituting  z=0  into  Eq.  (3.23),  we  can  obtain  the  temperature  at  the  end  of  the 
condenser  section  with  diffuse  effects  of  non-condensable  gases, 


4  =TS+ a  at  z  =  0  (3.29) 

Therefore,  the  temperature  drops  along  the  condenser  section  with  the  diffuse  effects  of 
non-condensable  gases  are 


A  Tn,=T-T,  = 


a[exp(mz)  - 1] 

Tl  -  {Ts  +  a  +  c[exp(m(2Lc  -  z))  -  exp(/wz)] } 


0<z<  Lcn 

^CJl  “  “  “ 


(3.30) 


Usually,  the  amount  of  non-condensable  gases  in  the  condenser  section  is  so  small 
that  their  diffuse  effects  can  not  reach  the  evaporator  section.  Accordingly,  the 
temperature  distribution  in  the  evaporator  section  is  not  affected  by  the  diffuse  of  non¬ 
condensable  gases,  and  Eq.  (3.7)  is  still  suitable  for  the  temperature  distribution  along  the 
evaporator  section.  But  the  temperature  distribution  in  the  condenser  section  will  be 
affected  strongly  by  the  diffusion  effects  of  non-condensable  gases.  At  the  same  time,  the 
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vapor  temperature  drop,  A Tr:,  caused  by  the  friction  and  centrifugal  force  must  be 

included  in  the  temperature  drop  in  the  condenser  section.  If  the  temperature  drop  caused 
by  the  diffuse  effects  of  non-condensable  gases  is  superposed  on  the  temperature  drop 
produced  by  vapor  friction  and  centrifugal  force,  the  following  relations  are  obtained 


&T:=ATr:+AT„: 

(3.31a) 

1 

E-h 

1 

ea 

II 

'ea 

1 

II 

•41 

E*"h 

<] 

(3.31b) 

T;„=Tv,L-ATv-ATnXc 

(3.31c) 

i 

n 

ii 

<3 

(3.31d) 

where  T,a  is  the  temperature  at  the  end  of  the  condenser  section  with  the  diffuse  effects 
of  non-condensable  gases  and  effects  of  the  friction  and  centrifugal  force  in  the  heat  pipe. 
Substituting  Eqs.  (3.31a),  (3.3 lb),  and  (3.3 lc)  into  (3.3 Id)  and  rearranging  it  yield 


Tv,  -  Tv  l  -  TLc  +  TZn  -  A Tr  +  A Tv  z  +  ATn: 


(3.32) 


It  should  be  pointed  out  that  the  temperature  drop  caused  by  the  vapor  friction  and 
centrifugal  force  must  be  dropped  from  Eq.  (3.32)  in  0  <  z  <  Lcn  since  the  vapor  does  not 
exist  in  this  section.  As  a  result,  after  the  vapor  temperature  drops  caused  by  vapor 
friction,  centrifugal  force,  and  diffuse  effects  of  non-condensable  gases  are  considered 
completely,  the  vapor  temperature  distributions  along  the  heat  pipe  length  are: 

In  the  evaporator  section. 


T  _T  TyFvQc  ( L-z )2  7>2sin^  1  1 

T„^-T„r  l  A#  l L  T  T  J}  7  [(Z0  +  ^  L)L  (z0  +  ^  z)z] 


L  v,L 


hjgPv  sin^ 


L-L, 


h 


fg 


<  z  <  L 


(3.33) 
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In  the  condenser  section, 


T„  -  TvL  - TLc  +Ts+a-  ATV  +  A Tr:  +  A Tn  ; 


{ ATv,z  -  ATv  +  Tv,l  -  c[exp(m(2Lc  -  z ))  -  exp(wz)]  Lcn  <  z  <  Lc 


(3.34) 


TrL  - TLc  +  Ts  +  a-  A Tr  +  n[exp(wz) - 1]  0  <z<Lc 


m  =  [ 


2  r„h 


£_£  12 


W-'i) 


(3.35) 


Eqs.  (3.33),  (3.34)  and  (3.35)  give  the  completely  closed-form  analytical  solution 
of  temperature  distributions  for  radially  rotating  miniature  heat  pipes  with  diffuse  effects 
of  non-condensable  gases.  A Tv  :  in  Eq.  (3.34)  is  the  vapor  temperature  drop  caused  by 
the  friction  and  centrifugal  force  along  the  direction  z,  which  is  given  by  Eq.  (3.4).  With  a 
low  and  moderate  rotation  speed,  ATr  .  is  very  small  and  negligible.  A Tr  in  Eq.  (3.34)  is 

the  total  vapor  temperature  drop  along  the  heat  pipe  length  without  the  diffuse  effects  of 
non-condensable  gases,  which  is  calculated  by  Eq.  (2.70).  For  the  heat  pipe  with  fins  in 
the  outer  surface  of  the  condenser  section,  r0  may  be  replaced  by  the  effective  radius  of 
outer  surface.  T  and  Ty  are  the  temperature  at  the  evaporator  end  cap  and  the  operating 
temperature  of  the  heat  pipe,  respectively.  In  the  normal  case,  T  can  be  replaced  by  T 
as  the  designed  heat  pipe  temperature  in  the  heat  pipe  design  and  calculation. 
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Vapor  Flow  Direction 


Fig.  3.1  Schematic  of  diffuse  effects  of  non-condensable  gases  in  the  condenser 
section. 


4.  HIGH-SPEED  ROTATING  TEST  APPARATUS  AND 


EXPERIMENTAL  PROCEDURES 

4.1  High-Speed  Rotating  Test  Apparatus  and  Radially  Rotating  Miniature  Heat 
Pipes 

In  order  to  prove  the  foregoing  analytical  results,  a  high-speed  test  apparatus  and 
a  data  acquisition  system  are  constructed.  The  data  acquisition  system  is  shown  in  Fig. 
4.1.  The  adjustable  revolution  range  of  the  rotating  test  apparatus  was  from  0  to  3,600 
rpm,  which  was  controlled  by  an  AC  inverter.  The  maximum  working  temperature  of  the 
miniature  heat  pipe  was  900  °C.  An  electric  heater  was  used  to  heat  the  heat  pipe.  Its 
electrical  power  was  adjusted  by  a  transformer  from  0  to  120  V,  and  supplied  by  a  two- 
channel  slip  ring.  The  temperature  distributions  along  the  heat  pipe  were  measured  by 
five  thermocouples,  and  fed  to  a  data  acquisition  system  through  a  5 -channel  slip  ring. 
All  experimental  data,  such  as  heat  inputs  and  temperatures  at  various  locations  of  the 
heat  pipe,  were  scanned  and  recorded  by  a  data  acquisition  system  (HP  34970A,  Data 
Acquisition/Switch  Unit),  and  the  data  were  saved  and  shown  in  a  computer 
subsequently. 

The  high-speed  rotating  test  apparatus  is  shown  in  Fig.  4.2.  For  this  particular 
apparatus,  the  main  parameters  which  must  be  controlled  in  the  experiments  were  the 
vibration  of  apparatus,  revolutions  of  the  test  apparatus,  temperatures,  heat  inputs,  and 
heat  dissipation  rate  of  the  condenser  section  along  the  radially  rotating  miniature  heat 
pipe.  The  shaft  was  driven  by  a  motor  [1]  with  1  horsepower  capacity,  in  which  the  motor 
revolution  was  varied  by  from  0  to  3,600  rpm.  In  order  to  dampen  the  vibration  of  the 
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rotor,  a  flexible  coupling  [18]  was  used  to  connect  the  shaft  with  the  motor.  Two  ball 
bearings  [6]  were  selected  to  fix  the  position  of  the  shaft.  The  partial  weight  of  the  rotor 
was  transferred  to  the  frame  [19]  through  the  bearing  house  [4]  and  the  4  sets  of  springs 
[3],  and  the  other  weight  was  transferred  to  the  frame  [19]  through  the  motor  [1]  directly. 
The  springs  [3]  connecting  the  rotor  and  frame  had  some  damping  function  when  the 
rotor  system  had  a  slight  unbalance  and  vibration.  Therefore,  the  rotating  test  apparatus 
had  some  self-balance  capability.  Two  slip-ring  assemblies  [7]  were  mounted  on  the 
shaft.  One  was  used  to  supply  the  electrical  current  from  the  power  transformer  to  the 
heater  [15]  that  heated  the  heat  pipe,  and  the  other  was  used  to  feed  the  temperature  data 
from  the  thermocouples  fixed  on  the  heat  pipe  to  the  data  acquisition  system.  The  rotor 
disk  consisted  of  an  inner  cylinder  [10]  and  an  outer  cylinder  [1 1],  The  inner  cylinder  was 
mounted  on  the  shaft,  and  the  outer  cylinder  was  fixed  on  the  inner  cylinders.  In  order  to 
prevent  the  outer  and  inner  cylinder  from  sliding  with  respect  to  each  other,  two  nuts  with 
a  bushing  were  bolted  between  the  inner  cylinder  and  the  outer  cylinder.  The  radially 
rotating  miniature  heat  pipe  [14]  and  a  counter  weight  [9]  were  mounted  on  the  opposite 
sides  of  the  outer  cylinder  in  order  to  balance  the  rotor  system.  The  heater  cap  [12] 
connected  on  the  outer  cylinder  had  two  functions:  one  was  to  protect  the  heater,  the 
other  was  to  change  the  cooling-air  flow  rate  in  the  condenser  through  the  change  of  the 
flow  window  of  cooling  air  at  the  heater  cap.  If  the  outer  diameter  of  outer  cylinder  was 
extended,  the  rotating  radius  of  the  heat  pipe  could  be  increased.  Thus,  the  high  speed 
rotating  test  apparatus  would  have  better  versatility  so  that  it  can  perform  a  range  of 
experiments  with  different  rotating  radii,  revolutions,  heat  inputs,  and  cooling  conditions 
for  the  radially  rotating  miniature  heat  pipe. 
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The  radially  rotating  miniature  heat  pipe  is  a  high-temperature  sodium-stainless 
steel  (Type  304W)  heat  pipe.  The  compatibility  of  304W  stainless  steel  shell  with  the 
sodium  working  fluid  has  been  well  documented  in  the  high-temperature  heat  pipe 
literature. 

Two  radially  rotating  miniature  heat  pipes  were  designed  to  operate  at  a  vapor 
temperature  of  about  800  °C  in  the  experiment.  One  was  1 .5  mm  in  inner  diameter,  and 
80  mm  in  length,  as  shown  in  Fig.  4.3.  The  other  was  2  mm  in  inner  diameter,  and  80  mm 
in  length,  as  shown  in  Fig.  4.4.  To  facilitate  the  sodium  charge  for  these  heat  pipes,  a 
reservoir  with  4  mm  in  diameter  and  3  mm  in  length  was  designed  on  top  of  the 
evaporator  section  with  an  inner  diameter  of  1 .5  mm.  The  other  reservoir  with  the  same 
diameter,  but  5  mm  in  length  was  designed  on  the  top  of  the  evaporator  section  with  an 
inner  diameter  of  2  mm.  The  end  cap  was  welded  on  the  reservoir  of  the  evaporator 
section,  and  a  filling  tube  with  an  outer  diameter  of  2.1  mm  was  welded  onto  the  end  cap. 
Two  radially  rotating  miniature  heat  pipes  had  the  same  outer  diameter,  6  mm,  and  the 
same  evaporator  and  condenser  lengths,  40  mm,  respectively.  Four  fins  were  made  on  the 
outer  surface  of  the  condenser  section  in  order  to  increase  the  convective  heat  transfer  on 
the  outer  surface  of  the  condenser  section.  The  heat  pipes  and  heater  caps  were  fabricated 
from  Type  304W  stainless  steel.  The  screw  thread  was  machined  on  the  heat  pipe  root  so 
that  the  heat  pipe  could  be  screwed  into  the  outer  cylinder  of  the  high  speed  rotating  test 
apparatus.  All  of  the  heat  pipe  parts  were  carefully  fitted  and  cleaned  according  to 
standard  procedures.  The  heat  pipes  and  end  caps  were  degreased  with  1,1,1- 
trichloroethane,  rinsed  in  tap  water,  then  rinsed  in  methanol  and  allowed  to  air  dry.  The 
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welding  process  of  filling  tubes,  end  caps  and  heat  pipes  were  protected  from  oxidation 
by  a  cover  layer  of  the  argon. 

The  sodium  charge  in  the  radially  rotating  miniature  heat  pipes  was  processed  by 
THERMACORE  INC.,  PA.  The  processing  procedure  was  as  follows: 

1)  Sodium  was  charged  into  the  heat  pipes  in  an  argon-filled  glove  box. 

2)  The  open  ends  of  the  heat  pipes  were  attached  to  a  turbo-pump  vacuum 
station. 

3)  While  pumped  at  a  high  vacuum  level,  heat  pipes  were  heated  to  100-200  °C 
to  drive  off  any  absorbed  water  and  atmospheric  gases. 

4)  The  filling  tubes  of  the  heat  pipes  were  pinched  and  welded  if  the  vacuum  in 
the  box  remained  at  a  high  vacuum  level. 

For  the  radially  rotating  miniature  heat  pipe  with  an  inner  diameter  of  1 .5  mm,  the 
sodium  charge  was  approximately  0.06  g,  and  0.08  g  of  sodium  was  charged  into  the 
other  radially  rotating  miniature  heat  pipe  with  an  inner  diameter  of  2  mm. 

There  were  five  thermocouples  of  type  K  used  to  measure  the  heat  pipe 
temperature  distributions  along  the  heat  pipe  length,  which  were  calibrated  with  a 
accuracy  of  +0.5  °C  by  the  OMEGA  ENGINEERING  INC.  Two  thermocouples  were 
installed  in  the  evaporator  section  to  measure  the  temperatures  in  the  evaporator  section, 
and  other  three  thermocouples  were  mounted  in  the  condenser  section  to  measure  the 
temperatures  of  the  condenser  section.  These  five  thermocouple  heads  were  inserted  into 
five  small  holes  of  diameter  of  0.8  mm  in  the  heat  pipe  shell  along  the  heat  pipe  length, 
and  were  fixed  firmly  on  the  outer  surface  of  the  heat  pipe. 
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The  heater  was  made  of  the  high  temperature  chemical  set  cement, 
OMEGABOND  “600”,  and  the  Nickel-Chromium  alloy  resistance  heating  wires,  NI80- 
010  or  NI80-012,  from  OMEGA  ENGINEERING  INC.  The  procedure  for  the  fabrication 
of  the  heater  was  as  follows: 

1)  Mixing  100  parts  of  the  high  temperature  chemical  set  cement  with  13  parts 
water  by  weight. 

2)  Painting  a  thin  layer  of  the  mixture  on  the  outer  surface  of  the  heat  pipe 
uniformly. 

3)  Wrapping  the  Nickel-Chromium  alloy  resistance  heating  wire  on  the  surface 
of  the  cement  as  a  heating  coil  after  the  high  temperature  chemical  set  cement 
cures  with  an  internal  chemical-setting  action  in  18-24  hours  at  ambient 
temperature  or  4  hours  at  1 80-220  °F. 

4)  Painting  another  thick  layer  of  the  mixture  on  the  heating  coil. 

The  heater  was  mounted  on  the  outer  surface  of  the  evaporator  section  firmly, 
and  the  insulation  material  of  about  16  mm  in  thickness  was  stuffed  between  the  heater 
and  heater  cap  to  make  sure  that  the  heat  from  the  heater  could  be  transferred  into  the 
evaporator  section  of  the  heat  pipe  effectively. 
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Heat  Pipe 


High-Speed  Rotating  Test  Apparatus 


Fig.  4.1  Schematic  of  a  high-speed  rotating  test  apparatus  and  data  acquisition 
system. 


Material:  304  Stainless  Steel 


Schematic  of  radially  rotating  miniature  heat  pipe  with  a  inner  diameter 


Material:  304  Stainless  Steel 


Schematic  of  radially  rotating  miniature  heat  pipe  with  a  inner  diameter  of  2  mm. 


4.2  Experimental  Results  and  Analyses 


According  to  the  foregoing  analysis  of  flow  and  heat  transfer  for  the  liquid  film 
and  vapor  in  a  radially  rotating  miniature  heat  pipe,  it  can  be  found  that  the  temperature 
drop  and  the  temperature  distribution  along  the  heat  pipe  length  are  the  function  of 
rotating  speed,  heat  pipe  length,  heat  pipe  inner  diameter,  working  fluid  properties,  tilt 
angle,  heat  input  in  the  evaporator  section,  and  cooling  condition  in  the  condenser 
section.  At  the  same  time,  other  parameters  such  as  non-condensable  gases  existing  in  the 
heat  pipe  and  heat  conductivity  of  the  heat  pipe  container  material  also  have  a  distinct 
influence  on  the  temperature  distribution  and  heat  transfer  in  the  heat  pipe. 

In  the  experiments  for  the  radially  rotating  miniature  heat  pipe,  the  tilt  angle,  <p , 
of  the  heat  pipe  is  designed  as  90 0 .  The  ranges  of  heat  inputs,  geometric  dimensions, 
flow  rates  of  cooling  air,  and  rotating  speeds  or  frequencies  are  as  follows: 

47  W< <2  <325  W,  d,= 1.5  mm  and  2  mm,  1=80  mm,  Lc=Le=A0  mm,  La  =  0, 

1.05x  10-3  m3/s<  W  <  13.4x  10"3  m3/s,  470  <  *>2Z0 /g  <  1881 ,  30 Hz  <  f<60Hz.\n 
the  following,  experimental  data  are  presented  in  terms  of  the  following  five  aspects. 
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4.2.1  Operating  Characteristics  of  High-Temperature  Radially  Rotating  Miniature 

Heat  Pipes 

As  mentioned  earlier,  there  are  many  influential  factors  for  the  operation  of  high- 
temperature  radially  rotating  miniature  heat  pipe.  The  major  factors  are  the  heat  input  in 
the  evaporator  section,  flow  rate  of  the  cooling  air  in  the  condenser  section,  heat  pipe 
size,  and  dimensionless  centrifugal  force  and  non-condensable  gases  in  the  heat  pipe. 

To  illustrate  the  effectiveness  of  the  heat  pipe.  Table  4.1  and  Fig.  4.5  show  the 
comparisons  of  the  temperature  distributions  along  a  heat  pipe  and  those  along  a  heat 
pipe  container  (the  heat  pipe  shell  without  charging  any  fluid)  with  the  same  flow  rate  of 
cooling  air,  dimensionless  centrifugal  forces,  and  geometrical  dimensions.  From  Table 

4.1  and  Fig.  4.5,  it  can  be  seen  that  the  heat  transfer  capacity  of  the  heat  pipe  container  is 
very  low,  which  is  47  W  for  the  heat  pipe  container  with  an  inner  diameter  of  2  mm,  and 
60  W  for  that  of  1.5  mm.  The  temperature  distributions  along  most  of  the  container 
length  are  approximately  linear,  and  temperature  in  the  portion  near  the  condenser  end  is 
close  to  the  temperature  of  the  cooling  air.  For  the  heat  pipes  with  the  same  flow  rate  of 
cooling  air,  dimensionless  centrifugal  force,  and  geometrical  dimensions,  much  more 
heat  can  be  transferred,  which  is  up  to  280  W  for  the  heat  pipe  with  an  inner  diameter  of 
2  mm,  and  250  W  for  that  of  1.5  mm.  If  the  thermal  conductivity  for  the  copper  is  386 
W/m'°C,  the  effective  thermal  conductance  of  the  heat  pipes  is  60  -  1 00  times  higher  than 
the  thermal  conductivity  of  copper.  These  comparisons  prove  that  high-temperature 
radially  rotating  miniature  heat  pipes  can  work  very  effectively.  Their  heat  transfer 
characteristics  are  much  better  than  any  metals. 
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Temperature  distributions  along  the  evaporator  section  and  most  portion  of  the 
condenser  section  are  almost  uniform.  It  shows  that  the  heat  pipe  works  perfectly  in  these 
sections.  But  there  exists  a  larger  temperature  gradient  near  the  end  cap  of  the  condenser. 
This  is  because  the  condenser  end  is  fixed  on  the  outer  cylinder  of  the  high-speed  rotating 
test  apparatus  with  a  screw  thread,  an  extra  heat  sink  is  formed  at  the  condenser  end,  and 
more  heat  from  the  condenser  end  is  transferred  into  the  heat  sink.  As  a  result,  because  of 
the  heat  sink  at  the  condenser  end,  the  temperature  distribution  close  to  the  condenser  end 
is  lowered. 

Table  4.2  and  Fig.  4.6  give  the  temperature  distributions  along  the  heat  pipe 
length  in  the  process  of  the  heat  pipe  startup.  At  r=0,  a  heat  input,  Q=176  W,  is  applied 
suddenly  to  the  evaporator  section  through  the  heater  described  early  for  the  heat  pipe 
with  an  inner  diameter  of  1.5  mm,  and  Q=175  W  for  the  other  heat  pipe  with  an  inner 
diameter  of  2  mm.  At  this  moment,  the  temperatures  along  the  heat  pipe  length  are  the 
same,  and  the  temperature  distribution  along  the  heat  pipe  length  is  a  straight  line  because 
of  the  retardation  of  the  heat  transfer  from  the  heater,  as  shown  in  Fig.  4.6.  With  the 
increase  in  time,  the  temperature  in  the  outlet  of  the  evaporator  section  rises  faster  than 
that  near  the  end  cap  of  the  evaporator,  and  temperatures  along  the  condenser  section 
almost  keep  constant.  Because  the  evaporator  end  has  a  larger  metal  mass  and  holds 
almost  all  the  solid  sodium,  it  needs  to  absorb  more  heat  to  raise  its  temperature.  In 
general,  this  phenomenon,  in  which  the  temperature  at  the  outlet  of  the  evaporator  section 
is  higher  than  that  at  the  end  cap  of  the  evaporator,  lasts  only  about  60  seconds.  After  that 
time,  the  temperature  at  the  end  cap  of  the  evaporator  is  higher  than  that  at  the  outlet  of 
the  evaporator  section,  and  temperature  along  the  condenser  section  increases  gradually 
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because  the  end  cap  of  the  evaporator  is  insulated  and  the  heat  from  the  outlet  of  the 
evaporator  section  is  transferred  to  the  condenser  section  by  the  axial  heat  conduction 
along  the  container  wall.  As  the  temperature  in  the  evaporator  section  is  further 
increased,  the  sodium  in  the  evaporator  section  is  evaporated,  vapor  flows  from  the 
evaporator  section  to  the  condenser  section,  and  the  heat  pipe  begins  working.  The 
temperatures  along  the  heat  pipe  length  rise  quickly.  After  about  1200  -1 800  seconds,  the 
temperature  distribution  along  the  heat  pipe  length  reaches  a  steady  state  and  the  process 
of  the  heat  pipe  startup  is  completed. 

At  a  lower  heat  input,  the  heat  pipe  could  not  work  since  the  operating 
temperature  in  the  evaporator  section  is  too  low  and  the  sodium  may  remain  to  be  solid 
and  would  not  evaporate.  When  the  heat  input  is  increased,  the  evaporator  section  of  the 
heat  pipe  starts  to  work  first.  But  the  vapor  generated  in  the  evaporator  section  will  be 
condensed  quickly  near  the  condenser  inlet,  and  condensate  will  return  to  the  evaporator 
section  under  the  action  of  the  centrifugal  force  if  the  flow  rate  of  the  cooling  air  in  the 
condenser  section  is  high  enough.  As  the  heat  input  in  the  evaporator  section  is  further 
increased,  the  front  of  vapor  flow  from  the  evaporator  section  will  move  forward  towards 
the  end  of  the  condenser  section,  and  the  working  section  in  the  condenser  will  be 
increased  gradually. 

Fig.  4.7  illustrates  the  steady  state  operation  of  two  radially  rotating  miniature 
high-temperature  heat  pipes  with  inner  diameters  of  1.5  and  2  mm  when  the 
dimensionless  centrifugal  force  is  470,  the  flow  rate  of  the  cooling  air  is  0.0067  m3/s,  and 
the  temperature  of  the  cooling  air  is  about  30  °C.  From  Fig.  4.7,  it  is  clear  that  the 
evaporator  sections  of  the  heat  pipes  start  to  work,  and  the  temperature  distributions  in 
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the  evaporator  section  are  almost  uniform  when  the  heat  input  reaches  100  W.  Because 
the  working  sections  in  the  condensers  are  relative  short,  the  temperatures  in  the 
condenser  section  drop  sharply  so  that  the  temperatures  near  the  condenser  end  are 
almost  close  to  the  cooling  air  temperature.  With  the  increase  of  heat  inputs,  the  operating 
temperatures  in  the  evaporator  sections  will  rise,  and  vapor  flow  rates  from  the 
evaporator  sections  will  increase.  At  the  same  time,  the  working  sections  in  the 
condensers  will  be  extended  towards  the  condenser  end,  and  the  temperature  along  the 
condenser  section  will  go  up  quickly.  When  the  heat  input  approaches  about  280  W  for 
the  heat  pipe  with  an  inner  diameter  of  2  mm,  and  250  W  for  that  of  1.5  mm,  the 
temperature  distributions  in  the  evaporator  sections  and  most  parts  of  the  condenser 
sections  are  almost  uniform.  However,  the  temperature  gradient  near  the  end  of  the 
condenser  section  still  remains  high  due  to  the  diffuse  effects  of  non-condensable  gases 
and  heat  sink  at  the  condenser  end. 
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Table  4.1  Comparisons  of  heat  pipes  and  heat  pipe  containers  with  the  same  geometry. 


(a)  co2Zal  g  =470,  dj=  2  mm,  W=6.7xl0'3m3  Is,  Q=280  W  for  the  heat  pipe 
and  Q=47  W  for  the  heat  pipe  container 


Dimensionless  heat  pipe 
length 

z/L 

Temp,  distributions  for 
the  heat  pipe 

rc) 

Temp,  distributions  for 
the  heat  pipe  container 

CC) 

0.000 

830.0 

814.2 

0.500 

821.8 

382.1 

0.750 

804.4 

79.5 

0.875 

742.4 

42.0 

1.000 

390.0 

37.9 

(b)  0)2Zal g=470,  <7,  =1.5  mm,  W=6.7xl0  ^m2  Is,  Q=250  W  for  the  heat  pipe 
and  Q=60  W  for  the  heat  pipe  container 


Dimensionless  heat  pipe 
length 

z/L 

Temp,  distributions  for 
the  heat  pipe 

CC) 

Temp,  distributions  for 
the  heat  pipe  container 

CC) 

0.000 

786.0 

801.5 

0.500 

776.4 

413.5 

0.750 

766.2 

91.8 

0.875 

698.4 

48.2 

1.000 

354.3 

35.3 
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(a)  co2Za!  g  =470,  d  =  2  mm,  W=6.7x  10  3  m3/s 


(b)  ru2Zfl/g=470,  d,  =1 .5  mm,  W=6.7 x  10~3  m3/s 


Fig.  4.5  Comparisons  of  temperature  distributions  between  the  heat  pipe 
and  the  heat  pipe  container. 


Table  4.2  Temperature  distributions  along  the  heat  pipe  length  with  the  diffuse 
effects  of  non-condensable  gases  during  the  heat  pipe  startup  process. 

(a)  co2ZJg  =  470,  d,  =1.5  mm,  0=176  W 


z/L=0 

z/L=0.5 

z/L=0.75 

z/L=0.875 

1! 

hJ 

T3 

t=0 

22  °C 

22  °C 

22  °C 

22  flC 

22  "C 

t=50  S 

28  "C 

54  °C 

22  °C 

22  "C 

22  "C 

t=60  S 

37  °C 

24  °C 

22  "C 

22  "C 

t=90  S 

329  "C 

244  °C 

37  °C 

25  "C 

22  "C 

t— 120  S 

379  °C 

61  °C 

32  °C 

23  °C 

t=150  S 

704  "C 

671  °C 

98  °C 

41  "C 

27  "C 

t=200 S 

794  °C 

792  "C 

564  °C 

211  °C 

81  "C 

t=l 800 S 

834  "C 

816  "C 

713  °C 

460  "C 

208  °C 

(b)  (o2Za/g  =470,  d=  2  mm,  0=175  W 


z/L=0 

z/L=0.75 

z/L=0.875 

z/L=l 

t=0 

23  "C 

23  °C 

23  "C 

23  °C 

23  "C 

t=50  S 

66  “C 

92  °C 

24  °C 

23  flC 

23  "C 

t=100  S 

| 

346  "C 

51  °C 

29  °C 

25  "C 

t=120 S 

438  °C 

70  "C 

35  "C 

26  "C 

t=140 S 

635  °C 

605  °C 

96  °C 

42  °C 

28  °C 

t— 160  S 

696  "C 

684  "C 

264  "C 

70  °C 

32  °C 

t=1200  S 

742  °C 

732  °C 

635  °C 

567  °C 

273  °C 
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-t=0 


-t=50S 
-t=60  S 
-t=90S 


H •—  t=120S 
— ♦— t=150S 
— ▲ — 1=200  S 
— — 1=1800  S 


(a)  co2Za/ g  =470,  4=1.5  mm,  0=176  W 


t=o  i 
t=50  S 
t=100S 
t=120S  | 
t=140  S  j 
t=160S  | 
t=1200  S | 


(b)  co2ZJg  =470,  d=2  mm,  0=175  W 

Fig.  4.6  Temperature  distributions  along  the  dimensionless  heat  pipe  length  with 
the  diffuse  effects  of  non-condensable  gases  during  the  heat  pipe  startup  process. 


Temperature  distributions  (  C)  Temperature  distributions  ( C) 


(a)  co2Za!  g  =470,  dj=2  mm,  W=6.7x  10  3  m3/s 


(b)  co2Za/g  =470,  J,=1.5  mm,  W=6.7x  10~3  m3/s 


Fig.  4.7  Temperature  distributions  along  the  dimensionless  heat  pipe  length 
with  different  heat  inputs. 


4.2.2  Effects  of  Heat  Inputs 


As  mentioned  earlier,  the  heat  input  has  a  very  important  influence  on  the 
temperature  distributions  and  heat  transfer  characteristics  of  the  heat  pipes.  When  the 
flow  rate  of  cooling  air,  dimensionless  centrifugal  force,  and  inner  diameters  of  the  heat 
pipes  are  fixed,  with  an  increase  in  heat  inputs,  the  operating  temperatures  in  the 
evaporator  section  will  rise,  and  the  temperature  distributions  along  the  condenser  section 
will  go  up  quickly,  as  shown  in  Figs.  4.8  and  4.9.  Accordingly,  the  heat  transfer 
characteristics  and  effective  thermal  conductance  will  be  improved  greatly.  If  the  flow 
rate  of  the  cooling  air  and  dimensionless  centrifugal  force  are  fixed  and  the  heat  input 
reaches  300  W  for  the  heat  pipe  with  an  inner  diameter  of  2  mm,  and  285  W  for  that  of 
1.5  mm,  the  temperature  distribution  along  the  heat  pipe  length  stays  relatively  uniform 
until  reaching  the  region  close  to  condenser  end.  But  a  temperature  peak  occurs  near  the 
condenser  end,  as  shown  in  Figs.  4.8(b)  and  4.9(a).  This  is  because  with  the  increase  of 
heat  inputs,  a  lot  of  vapor  is  generated  in  the  evaporator  sections.  If  the  flow  rate  of 
cooling  air,  dimensionless  centrifugal  force,  and  geometrical  dimensions  of  heat  pipes  are 
fixed,  more  vapor  would  flow  towards  the  end  of  the  condenser  section.  There,  the 
velocity  of  vapor  is  reduced  to  zero,  and  the  vapor  kinetic  energy  is  converted  into  the 
thermal  energy  which  causes  the  vapor  pressure  and  temperature  to  rise  again.  Owing  to 
the  effect  of  non-condensable  gases  and  the  heat  sink  at  the  condenser  end,  the 
temperature  at  the  condenser  end  is  still  relative  low.  As  a  result,  a  temperature  peak 
would  appear  close  to  the  end  of  the  condenser  section. 
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The  heat  transfer  rate  of  the  heat  pipe  can  be  changed  with  different  flow  rates  of 
cooling  air  and  dimensionless  centrifugal  forces  if  the  heat  input  and  geometrical 
dimensions  of  the  heat  pipe  are  kept  at  a  constant.  As  the  dimensionless  centrifugal  force 
and  flow  rate  of  the  cooling  air  are  increased,  the  film  thickness  of  condensate  on  the 
inner  wall  of  the  condenser  will  decrease  and  the  cross-section  area  of  the  vapor  flow 
passage  will  be  enlarged.  As  a  result,  the  heat  transfer  between  the  vapor  and  inner  wall 
of  the  condenser  will  be  enhanced.  At  the  same  time,  as  a  result  of  the  increased  flow 
rate  of  the  cooling  air,  much  more  heat  will  be  carried  away  from  the  condenser  section. 
Therefore,  the  vapor  from  the  evaporator  section  can  condense  more  effectively  in  the 
condenser  section  without  raising  the  temperature  in  that  region.  Figs.  4.8(c)  and  4.9(c) 
show  that  the  temperature  peak  for  heat  pipes  with  the  same  inner  diameters  has 
disappeared.  Compared  to  the  temperature  distributions  in  Fig.  4.8(b)  and  4.9(a),  the 
temperature  distributions  except  those  near  the  condenser  end  in  Fig.  4.8(c)  and  4.9(c) 
becomes  smoother  because  of  a  larger  flow  rate  of  the  cooling  air. 


(a)  o)2ZJg  =470,  d=2  mm,  W=6.7x  10'3  m3/s 


(b)  co2ZJg  =1306,  dj=  2  mm,  W=1 1.2x10  3  m3/s 


(c)  co2 Z a / g  =1881,  dt  =2  mm,  W=13.4x  10  3  m3/s 

Fig.  4.8  Temperature  distributions  along  dimensionless  heat  pipe  length  the 
with  different  rotating  frequencies  and  flow  rates  of  the  cooling  air. 
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(a)  co2Zal g  =470,  J(.=1.5  mm,  W=6.7xl0  3  m3/s 


(b)  (o2Zalg= 1306,  d,  =1.5  mm,  W=1 1.2x1 0‘3  m3/s 


„  900 
“  800 
I  700 
=  600 
1  500 

"S  400 

3  300 

03 

®  200 
E  100 
H  0 


(c)  o)2Zalg=\m,  4=1.5  mm,  W=13.4xl0'3  m3/s 

Fig.  4.9  Temperature  distributions  along  the  dimensionless  heat  pipe  length 
with  different  rotating  frequencies  and  flow  rates  of  the  cooling  air. 


0  0.2  0.4  0.6  0.8  1 

Dimensionless  heat  pipe  length,  z/L 
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4.2.3  Effects  of  the  Cooling  Air  Flow  Rates 


Flow  rate  of  the  cooling  air  in  the  condenser  section  has  a  dominant  influence  on 
the  operating  temperature  level  of  the  heat  pipe.  As  the  flow  rate  of  the  cooling  air  is 
increased,  the  convective  heat  transfer  on  the  outer  surface  of  the  condenser  section  will 
increased,  and  more  vapor  can  be  condensed  in  the  condenser  section.  Therefore,  the  heat 
transfer  capacity  of  the  heat  pipe  will  increase.  Accordingly,  if  the  heat  input  keeps  a 
constant,  with  the  cooling  air  flow  rate  in  increase,  the  operating  temperature  of  the  heat 
pipe  will  be  lower,  the  working  section  of  the  condenser  will  be  shorter,  and  the 
temperature  distributions  along  the  condenser  section  will  be  less  uniform  near  the 
condenser  end. 

Figs.  4.10  and  4.1 1  compare  the  temperature  distributions  along  the  dimensionless 
heat  pipe  length  with  the  same  heat  input,  dimensionless  centrifugal  force,  and  the 
cooling  air  temperature  when  the  flow  rate  of  the  cooling  air,  W,  is  changed  from  0.001 
m3/s  to  0.01 1  m3/s  through  the  change  of  the  flow  window  of  the  cooling  air  on  the  heat 
pipe  cap.  From  Figs.  4.10  and  4.11,  it  can  be  seen  that,  for  the  same  heat  input  and 
dimensionless  centrifugal  force,  the  temperature  distributions  along  dimensionless  heat 
pipe  length  rise  with  a  decrease  in  the  flow  rates  of  the  cooling  air.  Because  at  a  higher 
flow  rate  of  the  cooling  air,  the  operating  temperature  level  of  the  heat  pipe  is  reduced. 
Accordingly,  the  vapor  pressure  level  in  the  heat  pipe  is  reduced.  Due  to  this  reduced 
vapor  pressure,  the  viscosity  of  the  condensate  is  increased,  and  the  vapor  flow  resistance 
on  the  liquid-vapor  interface  will  rise.  Therefore,  the  temperature  drop  caused  by  the 
friction  will  be  increased.  At  the  same  time,  the  non-condensable  gases  in  the  heat  pipe  is 
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less  likely  to  be  compressed  towards  the  condenser  end  owing  to  the  reduced  vapor 
pressure.  As  a  result,  the  non-condensable  gases  would  occupy  a  larger  portion  of  the 
condenser  section,  the  temperature  distribution  near  the  condenser  end  becomes  less 
uniform,  and  the  temperature  gradient  is  steeper  near  the  condenser  end. 
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Dimensionless  heat  pipe  length,  z/L 


(a)  co2Za  /  g  =470,  d  =  2  mm,  Q=200  W 


(b)  o)2 Za  / g=  1306,  di=  2  mm,  Q=200  W 

Fig.  4.10  Comparisons  of  temperature  distributions  along  dimensionless 
heat  pipe  length  with  different  flow  rates  of  the  cooling  air. 
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(a)  co2  Za  I g  =470,  <=1.5  mm,  Q=200  W 


(b)  co2  Z a  I  g  =1306,  <  =1 .5  mm,  Q=200  W 

Fig.  4. 1 1  Comparisons  of  temperature  distributions  along  dimensionless 
heat  pipe  length  with  different  flow  rates  of  the  cooling  air. 
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4.2.4  Effects  of  Dimensionless  Centrifugal  Forces 


As  mentioned  in  the  foregoing  analyses  in  chapter  2,  at  a  higher  revolution  or 
rotating  frequency,  the  dimensionless  centrifugal  force  plays  a  dominant  role  in  the  total 
temperature  drop  along  the  heat  pipe  length.  With  an  increase  in  the  dimensionless 
centrifugal  force,  the  pressure  and  temperature  drops  along  the  heat  pipe  length  will  be 
increased  if  heat  inputs  and  geometrical  dimensions  of  the  heat  pipe  are  fixed.  Therefore, 
to  some  degree,  the  dimensionless  centrifugal  force  dominates  the  slope  of  temperature 
distribution  and  the  total  temperature  drop  along  dimensionless  heat  pipe  length. 

Though  the  dimensionless  centrifugal  force  has  a  similar  trend  for  the  total 
temperature  drop  along  the  heat  pipe  length  compared  to  the  flow  rate  of  the  cooling  air, 
the  mechanism  causing  the  change  of  the  temperature  drop  is  different.  As  the  cooling  air 
flow  rate  is  increased,  the  convection  on  the  outer  surface  of  the  condenser  section  will  be 
intensified,  the  working  section  in  the  condenser  will  be  shorted,  and  the  temperature 
distribution  along  the  condenser  length  will  become  steeper  if  the  heat  input  keeps  a 
constant.  But  the  slope  of  the  temperature  distribution  along  the  evaporator  section  will 
keep  a  constant.  As  a  result,  the  total  temperature  drop  along  the  heat  pipe  length  is 
increased.  While  the  dimensionless  centrifugal  force  is  increased,  the  temperature  drop 
caused  by  the  vapor  centrifugal  force  will  rise,  and  the  total  temperature  drop  along  the 
heat  pipe  length  will  be  increased.  Accordingly,  it  results  in  the  increase  of  the  slope  of 
the  temperature  distribution  along  the  heat  pipe  length.  Therefore,  the  analysis  here  is 
limited  to  the  slope  of  temperature  distributions  along  the  evaporator  section  in  order  to 
illustrate  the  effects  of  dimensionless  centrifugal  forces. 
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Figs.  4.12  and  4.13  show  temperature  distributions  along  dimensionless  heat  pipe 
length  when  heat  inputs  and  geometrical  dimensions  are  fixed  and  the  dimensionless 

centrifugal  force,  co2  Za  /  g ,  is  changed  from  470  to  1881.  As  shown  in  Figs.  4.12  and 
4.13,  at  a  low  rotating  frequency,  f=30  Hz,  which  corresponds  to  a  dimensionless 

centrifugal  force,  co2Za  I g ,  of  470,  the  slopes  of  the  temperature  distributions  along  the 
evaporator  section  are  relatively  small,  and  temperature  distributions  along  the  evaporator 
section  are  relatively  uniform.  When  the  rotating  frequency  and  dimensionless  centrifugal 
force  are  increased,  the  slopes  of  temperature  distributions  along  the  evaporator  section 
become  steeper  with  an  increased  temperature  drop  in  the  evaporator  section.  As  can  be 

seen,  at  a  higher  rotating  frequency  of  f-60  Hz,  corresponding  to  co2  Za  /g  =1881,  the 
slopes  of  the  temperature  distributions  and  the  temperature  drops  along  the  evaporator 
section  are  relatively  high.  With  a  further  increase  in  rotating  frequencies  and 
dimensionless  centrifugal  forces,  the  trend  will  be  more  evident. 
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(a)  Q=175  W,  di  =1.5  mm 


(b)  Q=200  W,  6?,.  =1.5  mm 


(c)  Q=225  W,  dj  =1.5  mm 

Fig.  4.13  Temperature  distributions  along  the  dimensionless  heat  pipe  length 
with  different  rotating  frequencies. 
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4.2.5  Effects  of  Inner  Diameters  of  the  Heat  Pipes 


The  inner  diameter  of  a  heat  pipe  determines  the  cross-sectional  area  of  vapor 
flow  in  the  heat  pipe.  For  a  radially  rotating  miniature  heat  pipe,  if  it  has  a  relatively 
larger  inner  diameter,  the  friction  of  vapor  flow  is  relatively  smaller.  As  a  result,  the 
temperature  drop  caused  by  the  friction  of  the  vapor  flow  is  smaller.  Conversely,  a 
radially  rotating  miniature  heat  pipe  having  a  smaller  inner  diameter  must  have  a  larger 
friction  of  vapor  flow  and  a  larger  temperature  drop  in  the  heat  pipe.  This  phenomenon 
has  been  analyzed  and  predicted  in  details  in  chapter  2.  The  experimental  evidence  will 
be  presented  in  this  chapter. 

Figs.  4.14  and  4.15  show  the  comparisons  of  temperature  distributions  along  the 
dimensionless  heat  pipe  length  for  heat  pipes  with  inner  diameters  of  2  mm  and  1.5  mm. 
It  is  evident  that  the  temperature  drops  along  the  heat  pipe  with  an  inner  diameter  of  1 .5 
mm  are  larger  than  those  with  an  inner  diameter  of  2  mm  when  the  heat  input, 
dimensionless  centrifugal  force,  and  flow  rate  of  cooling  air  are  kept  at  constants.  On  the 
other  hand,  the  operating  temperature  in  the  evaporator  section  of  the  heat  pipe  having  an 
inner  diameter  of  1.5  mm  is  higher  than  that  of  2  mm.  This  is  because  a  higher  pressure 
drop  is  needed  to  overcome  the  friction  of  vapor  flow  in  the  heat  pipe  with  an  inner 
diameter  of  1.5  mm.  Therefore,  the  vapor  pressure  in  the  evaporator  section  of  the  heat 
pipe  with  the  inner  diameter  of  1.5  mm  is  higher  than  that  of  2  mm.  Accordingly,  the 
saturated  temperature  in  the  evaporator  section  of  the  heat  pipe  having  an  inner  diameter 
of  1.5  mm  is  higher  than  that  of  2  mm.  This  trend  will  be  more  evident  with  a  further 
decrease  in  the  inner  diameter  of  the  heat  pipe. 
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(a)  Q=200  W,  co2Zal g  =470 


(b)  Q=200  W,  co2Zal g  =  1306 


(c)  Q=200  W,  a)2ZJg  =  m\ 


Fig.  4.15  Comparisons  of  temperature  distributions  along  the  dimensionless 
heat  pipe  length  with  the  same  flow  rates  of  the  cooling  air. 


4.3  Comparisons  of  the  Theoretical  Predictions  and  Experimental  Data  for  the 

Temperature  Distributions  along  the  Dimensionless  Heat  Pipe  Length 

According  to  the  analyses  in  chapter  3,  the  non-condensable  gases  have  an  important 
influence  on  the  temperature  distributions  along  the  heat  pipe  length.  The  analytical 
solution  based  on  Eq.  (3.7)  predicts  the  temperature  distributions  along  the  heat  pipe 
length  without  the  diffuse  effects  of  non-condensable  gases.  Eqs.  (3.33),  (3.34),  and 
(3.35)  give  the  closed-form  analytical  solution  of  temperature  distributions  along  the  heat 
pipe  length  with  the  diffuse  effects  of  non-condensable  gases.  These  analytical  solutions 
for  the  temperature  distributions  are  compared  with  the  experimental  results  under  the 
same  geometrical  dimensions  and  operating  conditions  of  the  heat  pipes. 

Table  4.2  and  Fig.  4.16  show  the  comparisons  among  the  analytical  solutions  without 
the  diffuse  effects  of  non-condensable  gases,  experimental  data  and  the  analytical 
solutions  with  the  diffuse  effects  of  non-condensable  gases.  It  is  clear  that  Eq.  (3.7)  gives 
a  set  of  ideal  analytical  predictions  without  diffuse  effects  of  non-condensable  gases.  In 
Eq.  (3.7),  the  temperature  distributions  along  the  dimensionless  heat  pipe  length  are 
caused  only  by  the  friction  on  the  liquid-vapor  interface  and  the  centrifugal  force  of  the 
vapor.  Therefore,  at  a  higher  operating  temperature,  the  temperature  drops  along  the 
dimensionless  heat  pipe  length  are  very  small,  and  the  temperature  distributions  are 
approximately  straight  lines,  as  shown  in  Fig  4.16.  The  experimental  results  under  the 
same  conditions  of  operating  conditions  and  geometrical  dimensions  are  used  in 
comparison  with  the  analytical  predictions.  It  can  be  seen  that  large  discrepancies 
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between  the  analytical  solutions  from  Eq.  (3.7)  and  experimental  results  are  present  in  the 
condenser  sections,  especially  close  to  the  condenser  ends. 

Eqs.  (3.33),  (3.34)  and  (3.35)  give  the  temperature  distributions  along  the 
dimensionless  heat  pipe  length  with  the  flat  front  of  vapor-gas  interface  when  the  length 
occupied  by  non-condensable  gases  in  the  condenser  section,  Lc  n ,  is  4  mm.  Compared  to 
the  closed-form  analytical. predictions  and  the  experimental  results  in  Table  4.2  and  Fig. 
4.16,  it  is  evident  that  the  analytical  predictions  from  Eqs.  (3.33),  (3.34)  and  (3.35)  are 
more  realistic  than  the  analytical  predictions  from  Eq.  (3.7)  and  are  more  close  to  the 
experimental  results.  According  to  Table  4.2  and  Fig.  4.16,  it  can  be  found  that  the 
analytical  solutions  from  Eqs.  (3.33),  (3.34)  and  (3.35)  agree  well  with  the  experimental 
results  in  the  entire  evaporator  section  and  most  part  of  the  condenser  section.  A  larger 
deviation  occurs  only  at  the  condenser  end.  This  is  because  the  condenser  end  cap  for  the 
tested  heat  pipe  forms  a  huge  heat  sink  since  it  is  connected  to  the  outer  cylinder  of  the 
rotating  experimental  apparatus,  and  most  of  the  heat  from  the  end  cap  of  the  condenser 
is  transferred  into  this  heat  sink  when  the  heat  pipe  operates.  As  a  result,  the  experimental 
data  from  the  end  cap  of  the  condenser  must  be  lower  than  the  predictions  from  the 
closed-form  analytical  solution. 

On  the  other  hand,  in  practical  cases,  for  the  radially  rotating  miniature  heat  pipe, 
the  vapor  condenses  on  the  inner  wall  of  heat  pipes,  which  forms  a  relatively  low- 
pressure  area  on  the  inner  wall  of  condenser  section.  The  non-condensable  gases  in  heat 
pipes  tend  to  move  toward  the  evaporator  section  under  the  action  of  the  centrifugal 
force.  Since  the  vapor  pressure  in  the  center  of  the  cross-section  is  higher  than  that  close 
to  the  inner  wall,  the  non-condensable  gases  will  diffuse  along  the  low-pressure  area 


93 


close  to  the  inner  wall.  As  a  result,  the  vapor-gas  interface  between  the  vapor  and  non¬ 
condensable  gases  is  not  a  flat  front,  but  a  spheroidal  one,  as  shown  in  Fig.  4.17.  This 
spheroidal  vapor-gas  interface  will  change  the  length  of  the  non-condensable  region,  and 
affect  the  heat  transfer  and  temperature  distributions  in  the  condenser  section.  Therefore, 
it  may  result  in  an  additional  deviation  between  the  closed-form  analytical  predictions 
and  experimental  data  in  the  condenser  section. 

Considering  the  complexities  of  the  flow,  heat  transfer  and  geometry  in  the 
condenser  section,  Eqs.  (3.33),  (3.34)  and  (3.35)  can  still  provide  a  good  prediction  for 
the  temperature  distributions  along  the  heat  pipe  length. 
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Table  4.3  Comparisons  of  analytical  solutions,  experimental  results  and  analytical 
solutions  with  the  diffuse  effects  of  non-condensable  gases  in  the  heat  pipe  (  Lc  n  =  4 
mm). 


(a)  co2Zalg  =470,  <=1.5  mm,  Q=250  W,  <=226  W/m2-°C 


Dimensionless  heat 
pipe  length 
z/L 

Analytical  solution  without 
non-condensable  gas 

("C) 

Analytical  solution  with 
non-condensable  gas 

rc) 

Experimental  results 

C’C) 

0.000 

786.0 

786.0 

786.0 

0.500 

782.1 

782.1 

776.4 

0.750 

779.1 

721.4 

766.2 

0.875 

778.5 

672.2 

698.4 

1.000 

778.2 

52i.8 

354.6 

(b)  (o2ZJg  =470,  <=2  mm,  Q=280  W,  <=297  W/m2-°C 


Dimensionless  heat 
pipe  length 
z/L 

Analytical  solution  without 
non-condensable  gas 

(°C) 

Analytical  solution  with 
non-condensable  gas 

rc) 

Experimental  results 

("C) 

0.000 

830.0 

830.0 

830.0 

0.500 

828.6 

828.6 

821.8 

0.750 

827.6 

764.3 

804.4 

0.875 

827.3 

711.9 

742.4 

1.000 

827.2 

551.4 

390.0 

(c)  co2ZJg=  836,  <=1.5  mm,  Q=280  W,  <=326  W/m2-°C 


Dimensionless  heat 
pipe  length 
z/L 

Analytical  solution  without 
non-condensable  gas 

CC) 

Analytical  solution  with 
non-condensable  gas 

CC) 

Experimental  results 

("C) 

0.000 

840.0 

840.0 

840.0 

0.500 

834.6 

834.0 

826.1 

0.750 

829.5 

771.5 

751.7 

0.875 

828.4 

718.0 

654.9 

1.000 

828.0 

549.7 

293.6 

(d)  co2Za/g=  836,  <=2  mm,  Q=280  W,  <=329  W/m2-°C 


Dimensionless  heat 
pipe  length 
z/L 

Analytical  solution  without 
non-condensable  gas 

CC) 

Analytical  solution  with 
non-condensable  gas 

CC) 

Experimental  results 

C'C) 

0.000 

819.0 

819.0 

819.0 

0.500 

816.4 

816.4 

802.0 

0.750 

813.6 

755.4 

752.5 

0.875 

813.1 

703.0 

639.0 

1.000 

812.9 

537.5 

314.6 
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(e)  co2Zalg=\m ,  J(.=1.5  mm,  Q=250  W,  4=336  W/m2-°C 


Dimensionless  heat 
pipe  length 
z/L 

Analytical  solution  without 
non-condensable  gas 

(°C) 

Analytical  solution  with 
non-condensable  gas 

re) 

Experimental  results 

re) 

0.000 

754.5 

754.5 

754.5 

0.500 

745.5 

745.5 

730.0 

0.750 

739.5 

680.3 

720.0 

0.875 

737.8 

642.2 

"520.0 

1.000 

737.2 

433.8 

235.7 

(f)  ofZalg  =  1881,  4=2  mm,  Q=280  W,  4=392  W/m2-°C 


Dimensionless  heat 
pipe  length 
z/L 

Analytical  solution  without 
non-condensable  gas 

(°C) 

Analytical  solution  with 
non-condensable  gas 

CC) 

Experimental  results 

("C) 

0.000 

680.0 

680.0 

680.0 

0.500 

665.3 

665.3 

661.2 

0.750 

654.4 

598.3 

608.0 

0.875 

651.6 

545.0 

577.5 

1.000 

650.7 

362.2 

248.9 
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Tempe 


Dimensionless  heat  pipe  length,  z/L 


Dimensionless  heat  pipe  length,  z/L 


(e)  co2Za  / g  =  1881,  <  =1.5  mm,  Q=250W 


(f)  co2Zalg  =  \m,  <=2  mm,  Q=275  W 

Fig.  4.16  Comparisons  of  analytical  solutions,  experimental  results,  and  analytical 
solutions  with  the  diffuse  effects  of  non-condensable  gases  ( Lc  n  =4  mm). 


Revolving  Shaft 


Condenser  Section  Vapor  Flow  Direction 


Fig.  4.17  Schematic  of  the  diffuse  effects  of  non-condensable  gases  at  a  higher 
revolution. 
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4.4  Error  Analyses  for  Experimental  Measurements  and  Analytical  Solutions 

Temperature  distributions  along  the  evaporator  length  exhibit  a  good  agreement 
between  the  analytical  solutions  and  experimental  data.  But  significant  deviations  exist  in 
the  condenser  section.  The  causes  for  these  deviations  are  very  complicated,  and  are 
relate  to  many  influential  factors.  An  error  analysis  is  provided  in  this  chapter,  and  the 
discrepancy  between  the  analytical  and  experimental  results  may  be  attributable  to  the 
following  factors: 

1.  The  errors  of  experimental  results,  which  include  the  uncertainties  in  the 
temperature  and  heat  input  measurements.  In  all  experiments,  the  accuracy  from  the 
thermocouple  of  type  K  is  0.5%,  the  accuracy  related  to  the  slip  ring  which  was  used  to 
transform  the  temperature  signals  is  about  0.5%,  and  the  HP34970  Data  Acquisition  has  a 
temperature  accuracy  of  about  0.2%.  In  addition,  since  it  was  difficult  to  weld  the 
thermocouples  onto  the  heat  pipe  shell,  a  high-temperature  chemical  cement  was  used  to 
fix  the  thermocouples  in  the  small  holes  drilled  in  the  heat  pipe  shell.  Therefore,  it  is 
possible  to  result  in  an  error  in  the  temperature  measurement  due  to  the  contact  thermal 
resistance  between  the  thermocouples  and  the  heat  pipe  shell  at  the  measuring  points.  If 
the  error  caused  by  the  contact  thermal  resistance  is  assumed  to  be  1%,  the  totally 
maximum  uncertainty  in  the  temperature  measurement  is  about  2.2%.  For  the  heat  input 
measurement,  the  heater  in  the  evaporator  section  is  insulated  with  an  insulation  material 
of  16  mm  in  thickness  stuffed  between  the  heater  and  the  heater  cap  in  order  to  reduce  the 
heat  losses  in  the  evaporator  section.  It  was  measured  in  experiments  that  the  temperature 
on  the  outer  surface  of  the  heater  cap  was  about  between  40-45  °C ,  which  corresponds  to 


101 


a  maximum  heat  loss  of  about  25  W  when  the  highest  temperature  in  the  evaporator 
section  was  between  800-850  °C  and  the  heat  input  was  250-300W.  On  the  other  hand, 
the  maximum  fluctuation  range  of  the  heat  input  in  experiments  was  +  2  W,  which  will 
result  in  an  error  of  about  0.8%  for  the  heat  input.  If  the  accuracy  of  the  wattmeter  and 
slip  ring  for  the  heat  input  is  about  1.5%,  the  total  uncertainty  in  the  heat  input  is  about 
1 0%.  As  a  result,  the  uncertainties  in  the  temperature  measurements  along  the  heat  pipe 
length  and  the  heat  input  in  the  evaporator  section  contribute  to  the  major  sources  of 
errors  in  experiments. 

2.  The  errors  of  analytical  solutions,  which  are  caused  by  one-dimensional  flow 
and  heat  transfer  modeling.  All  analytical  solutions  are  based  on  one-dimensional  flow 
and  heat  transfer  simplifications,  while  the  flow  and  heat  transfer  in  practical  heat  pipes 
should  have  three-dimensional  characteristics,  such  as  the  spheroidal  vapor-gas  interface 
at  the  condenser  section  at  a  high  rotating  speed.  Considering  the  miniature  flow  spaces 
and  diameters  in  radially  rotating  miniature  heat  pipes  and  complexities  of  the  flow  and 
heat  transfer,  one-dimensional  flow  and  heat  transfer  modeling  should  be  sufficient  to 
represent  the  major  characteristics  of  the  flow  and  heat  transfer.  However,  the  one¬ 
dimensional  simplification  will  cause  certain  deviations  between  the  analytical  and 
experimental  results. 

3.  The  errors  caused  by  experimental  apparatus,  which  include  the  interference 
produced  by  the  vibrations  of  experimental  apparatus  with  high  revolutions.  In  ideal 
cases,  it  is  reasonable  to  treat  the  flow  and  heat  transfer  problems  by  a  one-dimensional 
flow  and  heat  transfer  modeling  due  to  the  miniature  flow  spaces  in  radially  rotating 
miniature  heat  pipes.  But  in  practice,  as  the  revolutions  of  the  rotating  apparatus  was 
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increased,  relatively  large  vibrations  of  the  experimental  apparatus  were  created.  The 
one-dimensional  flow  and  heat  transfer  in  radially  rotating  miniature  heat  pipes  will  be 
interrupted  in  experiments  due  to  the  vibration  of  the  experimental  apparatus.  As  a  result, 
it  causes  three-dimensional  characteristics  of  the  flow  and  heat  transfer  in  experiments, 
which  enhances  the  deviations  between  the  analytical  solutions  and  experimental  results. 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  Conclusions 

The  radially  rotating  miniature  heat  pipe  is  an  effective  heat  transfer  device, 
which  uses  the  centrifugal  force  in  the  rotating  heat  pipe  to  return  the  condensate  from 
the  condenser  section  to  the  evaporator  section. 

In  this  report,  an  in-depth  study  for  the  radially  rotating  miniature  heat  pipe  has 
been  conducted  both  analytically  and  experimentally.  The  analytical  solutions  for  the 
condensate  film  and  vapor  flows,  film  thickness,  and  temperature  distributions  along  the 
heat  pipe  length  are  derived  based  on  appropriate  flow  and  heat  transfer  models.  The 
diffuse  effects  of  non-condensable  gas  on  the  temperature  distribution  along  the  heat  pipe 
length  are  also  studied.  Extensive  experimental  tests  on  radially  rotating  miniature  heat 
pipes  with  different  influential  parameters  are  conducted,  and  various  effects  on  the 
operation  of  the  heat  pipes  based  on  these  influential  parameters  are  searched.  Analytical 
solutions  for  the  temperature  distributions  along  the  heat  pipe  length  with  diffuse  effects 
of  non-condensable  gases  are  obtained  based  on  extensive  experimental  results. 
According  to  the  analytical  and  experimental  studies  presented  in  this  report,  the  follow 
conclusions  may  be  made: 

1.  The  radially  rotating  miniature  heat  pipe  has  a  simple  structure,  very  low 
manufacturing  cost,  and  very  high  heat  transfer  capability  and  effective  thermal 
conductance  (60-100  times  higher  than  that  of  the  copper).  At  the  same  time,  the  heat 
pipe  can  withstand  strong  vibrations,  and  work  in  the  high-temperature  and  adverse 
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environments.  Therefore,  the  combination  of  traditional  air-cooling  technology  with 
radially  rotating  miniature  heat  pipes  is  a  feasible  and  effective  cooling  means  for  high- 
temperature  turboengine  blades. 

2.  At  a  relatively  high  rotating  speed  or  a  large  revolving  radius,  co2  Z  /  g  >  150, 
the  maximum  temperature  drop  across  the  liquid  film  is  negligible.  However,  a  smaller 
heat  pipe  diameter  or  a  lower  heat  pipe  operating  temperature  may  increase  the 
temperature  drop.  Still,  in  any  cases,  the  maximum  temperature  drop  across  the  liquid 
film  is  much  less  than  1  K.  It  means  that  the  commonly  encountered  condenser  limitation 
for  low  temperature  heat  pipes  would  never  be  a  serious  problem  for  the  high  temperature 
heat  pipes. 

3.  A  heat  pipe  with  sodium  as  the  working  fluid  would  have  a  smaller  vapor 
temperature  drop  compared  to  the  heat  pipe  with  potassium  as  the  working  fluid.  For  the 
case  without  diffuse  effects  of  non-condensable  gases,  the  vapor  temperature  drop  is  on 
the  order  of  1 0-40  K.  This  temperature  drop  is  relatively  high  for  a  low-temperature  heat 
pipe.  However,  for  the  high-temperature  heat  pipes  studied  in  this  report,  this  temperature 
drop  is  small  compared  to  the  average  heat  pipe  operating  temperature  level. 

4.  The  inner  diameter  of  the  heat  pipe  determines  the  cross-section  area  of  the 
vapor  flow  in  the  heat  pipe.  For  a  radially  rotating  miniature  heat  pipe  with  a  larger  inner 
diameter,  the  friction  of  the  vapor  flow  is  smaller.  As  a  result,  the  temperature  drop 
caused  by  the  friction  of  the  vapor  flow  is  smaller.  Conversely,  a  radially  rotating 
miniature  heat  pipe  with  a  smaller  inner  diameter  would  have  a  larger  friction  of  the 
vapor  flow  that  would  cause  a  larger  temperature  drop  along  the  heat  pipe  length. 
However,  for  a  heat  pipe  with  dJL  greater  than  0.02  or  d,  greater  than  1.5  mm,  the 
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influence  of  the  heat  pipe  diameter  on  the  operating  temperature  level  and  the  vapor 
temperature  drop  is  relatively  small. 

5.  The  heat  input  has  a  very  important  influence  on  the  temperature  distributions 
and  heat  transfer  characteristics  of  a  heat  pipe.  When  the  flow  rate  of  the  cooling  air, 
dimensionless  centrifugal  force,  and  inner  diameters  of  the  heat  pipe  are  fixed,  with  an 
increase  in  heat  input,  the  operating  temperatures  in  the  evaporator  section  will  rise,  the 
working  section  in  the  condenser  will  be  increased,  and  the  temperature  level  in  the 
condenser  section  will  go  up.  Accordingly,  the  heat  transfer  characteristics  and  effective 
thermal  conductance  will  be  improved  greatly.  However,  a  radially  rotating  miniature 
heat  pipe  with  too  high  heat  input  will  encounter  its  heat  transfer  limitation. 

6.  The  dimensionless  centrifugal  force  plays  an  important  role  in  the  temperature 
drop  along  the  heat  pipe.  As  the  rotating  speed  of  a  heat  pipe  is  increased,  the 
dimensionless  centrifugal  force  will  rise,  and  the  pressure  and  temperature  drops  along 
the  heat  pipes  will  increase  accordingly.  At  a  low  rotating  speed,  the  contribution  of  the 
centrifugal  forces  is  relatively  small,  but  at  a  high  rotating  speed,  the  contribution  of  the 
centrifugal  forces  becomes  a  dominant  factor.  This  is  especially  obvious  for  a  heat  pipe 
having  a  relatively  large  diameter.  However,  at  a  fairly  high  rotating  speed  (15,000  rpm- 
20,000  rpm),  the  vapor  temperature  drop  caused  by  the  centrifugal  force  is  still  relatively 
small. 

7.  The  flow  rate  of  the  cooling  air  in  the  condenser  section  has  a  dominant 
influence  on  the  operating  temperature  of  the  heat  pipe.  When  the  flow  rate  of  the  cooling 
air  is  increased,  the  heat  transfer  capacity  of  the  heat  pipe  will  be  improved,  the  operating 
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temperature  of  the  heat  pipe  will  be  reduced,  the  working  section  of  the  condenser  will  be 
shorten,  and  the  temperature  distributions  near  the  condenser  end  will  become  steeper. 

8.  The  non-condensable  gases  in  the  heat  pipe  have  an  adverse  effect  on  the 
operation  of  the  heat  pipe.  The  non-condensable  gases  will  be  driven  to  the  condenser 
end  by  the  vapor  from  the  evaporator,  and  prevent  the  vapor  to  fill  the  entire  condenser 
section.  As  a  result,  the  non-condensable  gases  in  the  heat  pipe  will  result  in  a  large 
temperature  drop  near  the  condenser  end  and  reduce  the  thermal  conductance  and  heat 
transfer  performance  of  the  heat  pipe. 
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5.2  Recommendations 


1 .  The  radially  rotating  miniature  heat  pipes  can  be  made  with  different  geometrical 
cross-sectional  areas,  such  as  square,  triangular  and  elliptical  cross-sectional 
areas  in  order  to  meet  the  requirement  for  various  applications  in  industries. 

2.  Improving  the  performance  of  the  rotating  test  apparatus,  and  conducting 
experiments  for  radially  rotating  miniature  heat  pipes  at  a  higher  rotating  range 
that  is  close  to  the  operating  revolution  of  a  gas  turbine. 

3.  Seeking  a  better  high-temperature  chemical  cement  for  the  heater  so  that  higher 
heat  inputs  can  be  provided  in  experiments. 

4.  Improving  the  filling  procedure  for  radially  rotating  miniature  heat  pipes  so  that 
non-condensable  gases  can  be  eliminated  in  the  heat  pipe. 
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